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ABSTRAC T 

The study of the intervaiing fluid in the separation 
gap between bearing surfaces in supporting the imparted load 
has been a fundamental consideration in the field of lubrication. 
The quality, behaviour, range of utility in the physical 
environment together with the molecular structure of the 
lubricant and its affinity to the bounding surfaces important 
in understanding the lubricating mechanism in the perfozrmance 
of modern machinery. 

The need for fairly good theoretical prediction of bearing 
characteristics is possible when the realistic lubricating 
conditions such as rheological anomalies or variation of 
viscosity/ consistency in the fluid film., surface roughness, 
elastic deformation under heavy loads etc. are included in the 
analysis . 

It appears that the properties of the lubricant film in 
contact with the solid surface are different from that of the 
bulk lubricant' Kingsbury [ 1 ] was one of the earliest investigat- 
ors to report the rheological abnormalities of the fluid in 
the proximity of the solid surface- He observed an enhanced 
viscosity in the region of attraction of surface molecules of 
the. metals. The most important and interesting work in the 
study of rheological abnormalities of the fluid very near to 
the bearing surface is the work of Needs [2 ] who squeezed 
fluids between two steel plates and observed a stable residual 
film. More recently, the works of Cameron and his Cbworkers [3,4] 
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indicated an increase in the viscosity with the addition of 
certain additives and this results in a much larger viscosity 
than the bulk of the fluid. The presence of additives in the 
base oil characterises a non-Nevjtonian behaviour. The non-N evrtoni ^ 
behaviour of the lubricant raay be characterised by the power 
law model. 

Another important factor to affect the bearing performance 
is the ±.he3rmal affect. This effect is generally studied by 
coupling Reynolds and energy equations. The complexity of 
the complete solution to Reynolds equation arises mainly due 
to the inclusion of energy equation in the calculation procedure- 
A different approach to avoid energy equation was presented 
by Tipei and coworkers [5^.6] who replaced the viscosity-temperatu; 
relation by a film thickness-viscosity relation. With this 
approach / an extensive parametric study of lubrication problems 
especially in the case of power law lubricants may be attempted. 
This approach may be used to study the consistency variation ; 

along the film thickness. 

Prasad £ 7 ] analysed the lubrication of rollers using . a ; 

power law fluid considering consistency variation across the 
film thickness- Soma of his results obtained are open to question ; 
as the load ratios defined to study consistency variation across 
the film thickness was considered by taking into account cavitatiq 
and no cavitation load . In the present work a more general I 

model to account for consistency variation across.and along the fj 
thickness is developed to study the effect of this variation on tl 
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lubrication characteristics. Consistency variation across 
the film thickness is accounted through a consistency ratio 
factor while its variation along the film thicknessis accounted 
through a thermal factor- An attempt is made in the present 
v/ork to use this model for consistency variation in lubricating 
conditions-.’ The model is used for smooth as well as rough 
bearings and rigid as well as deforming bearings under heavy 
loads. The thesis work is divided into seven Chapters. 

Chapter 1 is introductory in nature. It contains a selected 
list of experimental and theoretical results in support of 
rheological anomalies in the convicinity of the bearing surface- 
Apart from the natural phenomenon of this rheological 
behaviour, a brief history of works relating to surface roughnes 
non-hev/tonian behaviour of the lubric-nts, elastohydrodynamics 
lubrication are dealt with in conformity with the theoretical 
analysis. 

In Chapter 2 , the generalized Reynolds equation considering 
consistency variation across as well as along the film 
thickness is derived. Subsequently, it has been used to study 
the effect of consistency variation in the case of rigid rollers 
Using Grubin^ s theory [ 8 ] the model for consistency variation 
is used in the inlet analysis of elastohydrodynamic lubrication. 
It has been established that the effect of the high consistency 
peripheral layer is to increase the load capacity in the case 
of rigid rollers and to increase the minimum film thickness 
in the ehd lubrication. In the later case , it was theoreticall 



established that no set pattern of minimum film thickness 
V 7 ith respect to the flov; behaviour index is observed^ When 
the total peripheral layer thictoess is of the same order 
as that of the minimum nominal film thickness. This is 
atributed to the complexity of the rheological effects in 
the ehd lubrication. The effect of the thermal factor in the 
case of high consistency peripheral layer is to decrease the 
load capacity in the case of rigid rollers and to decrease 
the minimum film thickness in the ehd case. 

The model for consistency variation is used in the analysis 
of squeeze film in Chapter 3. It is observed ttiat the presence 
of high consistency peripheral layer increases the load capacity 
and response time. The increase is more for pseudoplastics 
as compared to dilatants. It was concluded that consistency 
in the peripheral region and surface asperities, it was found 
that there exists a critical value for the consistency ratio, 
above which the trend of the load ratios with respect to the 
flow behaviour index are similar to that for the smooth case. 
When the consistency is higher in the central layer as ; ' 
compared to that in the peripheral layer, there is an increase 
in the load capacity for longitudinal roughness. For a high 
consistency peripheral layer, in the transverse roughness there 
is greater load capacity and protection against seizure. 

In Chapter 6 is studied the squeeze films on rough bearings 
considering consistency variation.. It is concluded that in the 
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presence of surface asperities, the effect of consistency 
variation across the film thickness is similar to the case 
of smooth bearing- The effect of roughness on load and 
response time is more for dilatants than pseudoplastics. 

In Chapter 7 is studied the effect of consistency variation 
on the inlet pressure for both longitudinal and transverse 
roughness. It is inferred that an increase in the consistency 
in the peripheral layer results in an increase in the inlet 
nominal film thickness. In the case of transverse roughness 
the inlet pressure increases as the value of the roughness 
parameter increases for both the cases v/hen the total peripheral 
layer thickness on both the rolling surfaces is greater or less 
than the minimum nominal film thickness. The effect of the 
thermal factor is to decrease the inlet pressure for transverse 
case and to increase it for longitudinal case. 
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CHAPTER 1 


GENERAL INTRODUCTION 
1.1 EARLY DAYS OF LUBRICATION ’ 

The manner in which the intervening fluid in the 
separation gap between the bearing surfaces supports 
the imparted load has been the subject of intensive study 
to understand the lubrication mechanism between the various 
moving parts of the modern machinery from gigantic 
structures to minute components of mechanical system. The 
vital study associated with the lubrication process centres 
around friction and wear of the machinery during running - 
in. The inevitable occurance of these natural phenomena cause 
severe operational disturbances and in some cases lead to 
total destruction of the entire mechanical system. One 
of the persistent demands, naturally, is to control friction 
and wear and mitigate wastage of energy during operation 
of the systona. To meet this demand, the study of the quality, 
behaviour, range of utility in the physical environment 
of the lubricating fluid is essential- Further, the molecular 
structure of the Ixibricant and its affinity to the 
bounding surfaces together with the strctural and topographical 
changes of these surfaces .. makes the study useful for 
practical applications. 

The lubrication practice is old. During the era of 
Egyptian Pharaohs, it was found that the axle of the chariot 
wheel had beoa lubricated with a sticky and slightly greasy 
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substance [. 1 ] • '^e substance contained road dirt such 

as quartz sand and ccanponents of aluminium and lime* Its 
melting point was determined to be 20‘^F. It was suggested 
that it could be either mutton or beef tallow, both of 
which were suitable for axle lubrication in that tropical 
country. Until very recently the lubricants have been of 
animal, vegetable or marine origin - mutton tallow, lard, 
goose grease, castor, cotton seed and other vegetable oils. 
The use of mineral oils dates back to a little over 
hundred years. Today straight mineral oils cannot cope 
up with the situations when the growing danand to transmit 
greater power through mechanisms smaller in size and weight 
imposes increasing burden on the operating lubricant. 
Recently, the synthetic lubricants, with greater thermal 
and viscosity properties than those of the mineral oils, 
are used especially in most arduous and transient bearing 
performances . 

Pioneering work in the field of tribology or lubrication 
was initiated by Adams [ 2 ] in 1853, while developing 
certain patents for railway axle bearings i However, the 
new science of lubrication on hydrodynamic principles was 
founded by Petrov [. 3 ] who theoretically analysed the 
friction effects and film lubrication. In the same year 
Petrov published his theory. Tower [ 4 ] dononstrated that, 
a loaded oil lubricated bearing might be subjected to local 
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pressure substantially higher than its mean pressure* Three 
years later, Reynolds £ 5 ] independent of Petrov's theory, 
was able to confirm Tower's results using the properties 
of the laminar fluid film in the equation which was named 
after him. Certain solutions to this equation were obtained 
by him that agreed well with Tower's experimental measurOTents. 
Later on, Sommerf eld 1 6] in 1904 made an elegant theoretical 
extension to journal bearing lubrication problems* 

To obtain solutions to Reynolds equation in lubrication 
probl^s , assumptions such as no slip, constancy of viscosity 
of the lubricant across and along the film thickness etc., 
are usually made. However, for bearings operating under 
actual environmental and physical conditions these assumptions 
may not be valid and hence, the prediction of bearing 
characteristics based on these theoretical -studies is'bound to 
differ with experimental measurements and observations. To 
predict the bearing characteristics with a fair degree of 
accuracy close to the actual performance, it is, therefore, 
indispensable to make generalisations of the Reynolds equation. 
Accordingly, various generalisations taking into account 
of the viscosity variation across and along the film thickness, 
thermal effects, porosity and slip condition, non-Newtonian 
characteristics of the lubricant, surface roughness, rheological 
anomalies of the lubricant very near to bearing sur.faces, 
elastic deformation of the bearing suefaces due to a highly 
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developed pressure etc. are made. Keeping these in view, 

lubrication characteristics of various designs of hybrid 

two lobe and three lobe bearings, multilayed porous bearings, etc> 

apart from the regular geometrical conventional bearings, 

have been determined. A historical and detailed review of 

the same would be an Herculean task and beyond the scope 

of the present work. 

1. 2 ■ LUBRICATION' AND RHEOLOGY 

It has been pointed out earlier that the properties of 
the lubricant film contiguous to the bearing surface are import- 
ant -in determining v the bearing characteristics. It 
appears that the properties of the liquid in contact with 
a mating solid surface, influenced by the proximity of the 
liquid to the bearing surface are a subject of intensive 
study and a matter of controversy. The molecules of the 
liquid near to adhering surfaces or of the adsorbed layer 
on the surface are believed to behave differently from the 
adjacent molecules of the liquid in bulk. The nature of 
the surface influence and depth of its penetration into the 
surface led to a variety of divergent opinions* It is generally 
accepted that the surface zone is not merely a monomolecular 
layer but a multimolecular layer i.e., the molecular orientation 
extends effectively to many molecular lengths £7,8]. Adamson 
[ 9 ] pointed out that though evidence of a deep surface 
orientation was more or less circxrastantial, the molecular 
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orientation does occur probably as a result of perturbations 
passed on from molecules to molecules. 

Several experiments were conducted to examine whether 
the fluid in the vicinity of the adhering surface developed 
properties significantly different from that of the fluid 
in bulk. Kingsbury 10 ] was one of the earliest investigat- 
ors to report the rheological abnormalities of the fluid in 
the proximity of the sold surface. He observed an enhanced 
viscosity of the fluid in the region of attraction of the 
surface molecules of the metals. This enhancement might 

range from a small increase in viscosity to rigidity, and the 

-7 -3 

depth of the affected zone might range frcxn 10 -10 mm. 

The most important and interesting work in the study of 
rheological abnormalities of the fluid very near to the 
bearing surfaces is the work of Needs [ 11 ] who squeezed 
fluids between two steel plates and observed a stable residual 
film.t Though some of his observations are questioned, 
it appears reasonable that an alteration of liquid film occurs 
near the bounding solid surface which makes it substantially 
more viscous compared to the bulk of the liquid. 

The most telling support to Needs observation is, perhaps, 
provided by the findings . of Fuks [12,13] who conducted a 
number of experiments using an apparatus similar to that used 
by Needs. It is of interest to note that an addition of 
0.1 ■/. of stearic acid to a chromatographic ally separated’ 
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naphthene-paraffin fraction oil has increased. t±ie thickness of 
the residual film. It is equally int^esting to note that 
lower Tnolecular weight hydrocarbons such as benzene, hexane, 
cyclohexane, isooctane and decane exhibited no residual 
'film* It was observed that the film thickness increased 
with the fatty acid concentration,, with surface energy of the 
solid substrate and with increasing hydrocarbon chain length 
of both solvent and fatty acid. 

The works of Needs and Fuks were a matter of great controversy 
and criticism by many workers, notably Hayward and Isdale £ 14 ] . 
It has been suggested that the asperities on the solid surface 
m.ight account for residual film. The dirt and the unsolved 
material in the fluid film have been also offered as an 
explanation. Drauglis et. al> £ 15 ] observed that if the data 
of Fuks were to be explained on the basis of dirtv then we 
had to postulate that that dirt would change its properties 
with temperature, length of additive molecule and concentration 
of additives, but such a dirt would be a most xinusual material. 
Derjaguin [16, 17] pointed out that the molecular orientation 
'' in depths would occur near solid boundaries and give rise to 
significant rheological changes- Increase in viscosity was 
observed where traces of polar or surface active substances 
were present. 

More recently, Cameron and his coworkers [ 18-21] studied 
the influence of surface active components in various lubrication 
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systOTS.The lubricants used were paraffins m=5inly hexadecane, and 

the additives consisted of various long chain polymer compoxinds. 
It was pointed out that the oil film was much influenced 
by the presence of additives. The scuffing load was found 
to increase with the chain length of additive and solvent. 
Squeeze film studies similar to Fulcs were made. It was 
observed that the plates came to rest at a separation 
of 2xl0~ mm. with all fluids including pure cetane- 
The increase in surface viscosity was greatest when additives 
were matched with carrier oils in terms of chain ler^th 

and shape. The results of these studies are explicable 

-5 —3 

by postulating the presence of a surface film 10 -10 mm. 

thick which has a much larger viscosity than the bulk of 
the fluid. 

Prom these studies we conclude , therefore# that the 
evidence of the existence of rheological abnormalities (like 
enhanc orient of viscosity) in the neighbourhood of solid 
boundaries on the whole seems to be positive, despite some 
controversy. It is clear that zones of enhanced viscosity 
(giving rise to residual film) would have important effects 
on the lubrication of practical machinery. 

As has been pointed out earlier, the presence of additives 
in base oils gives rise to rheological abnormalities in thin 
films such as enhancoaent of viscosity. Further, the 
incorporation of additives to the base lubricant improves; the 

/ 
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primary functions of the lubricant in minimizing sliding 
friction and offering surface protection against damage 
and wear [ 22 ]- Ordinary mineral lubricating oils contain 
various polar com pounds, many of them are surface active 
in a non-polar hydrocarbon mediim ,and thus belong to the 
class of liquids for which the eviderK:e of the effect of 
viscosity enhancement is strongest. The rheological anomalies 
could /therefore, have a considerable effect on the minimum 
film thicloiess/ which in turn influences wear and fatigue of 
machine components. 

The idea of blending additives to petroleum oils to 
improve the performance was used in the power, indxjstrial 
and railway lubricants. Modern lubricant additives -are based 
on many years of scientific research, designed to meet 
the extreme demands of modern machines and for high performance 
ratings under actual working conditions. The automative 
industry and machine equipment manufacturers use additives 
fortified lubricants to meet severe machine service dsnands, 
thereby widening the margin of safety of operations. The 
additives are incorporated in carrier oils for a variety of 
purposes. There are additives used as rust inhibitors 
(Amine phosphates), corrosion inhibitor (sulphurized olefins), 
fire resistors (hologenated hydrocarbons), detergents (calcium/ 
barium sulphorates) viscosity improvers (polymethylacrylate, 
powders of graphite and molybdenum disulphide), extreme 
pressure (ep) additives (sulphurised fats) etc. There are 
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certain experimental reports to support that the addition 
of additives to base oils reduce friction^ wear and oihance 
the lubrication process [ 23 ] . In particular, it has been 
noted that in the lubrication of rollers, the film thickness 
considerably increases when additives are introduced in the 
base oils [ 24 ]• The additives or siirf actants, thus, present 
in the lubricant enhances the effective viscosity of the 
lubricant near the contiguous siirface which is beneficial 
to the hydrodynamic Ixibrication process | 25 ] . 

The use of EP additives in lubrication is important-. It 
is usually accepted that these additives react with sliding 
metal surfaces to form a surface film which reduce wear and 
aids seizure resistance. The reaction between EP additives 
and sliding surfaces is promoted mainly by the increase in 
the surface temperature owing to frictional heat [ 26 ]. 

For extreme t^perature applications, synthetic liquids such 
as silicones, esters are used to improve temperature properties 
The thickening agents are usually lithium, calcium, sodium, 
barium and aluminium for temperature stability £ 27 ] . A 
more detailed knowledge than the existing one at present 
about the mechanism of formation and chemical structure 
of the reaction layers is of great interest for the development 
of new and more effective additives £ 28 ] . 
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1.3 SURFACE ROUGHNESS 

The relation between surface finish and wear is of 
great practical significance to bearing design engineers. 

The wear bdiaviour of the machine surface is affected by 
surface topography and surface integrity. The surface 
roughness and the surface waviness are two topographical 
structures, respectively, describing short range and long 
range geometric deviations of a surface from the nominal 
geometric shape [ 29 ]. Several attempts have been 
made to study the effects of roughness following deterministic 
and stochastic approaches. In deterministic approach, 
the film thickness is nornally assumed to be xepras^tod 
by sine ( cosine ) function for a given roughness profile 
and usual hydrodynamic equation is, then, solved to 
calculate the bearing characteristics [ 30 ] . It was 
showed that load capacity, frictional force etc. are 
dependent on the amplitude and wave length of these 
waves representing the roughness profile [ 31 ]- In 
the stochastic approach, the film thickness is assumed 
to be a stochastic or ergodic function and stochastic 
average is taken either for the Reynolds equation or 
for its solution. Chen and Sun [ 32 ] pointed out that 
the method of averaging was practical only for the 
analytical solution of the Reynolds equation . Tzeng and 
Saibel [ 33,34 ] introduced the stochastic concept of 
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surface roughness and obtained mean pressures for infinite 
slider and short joxirnal bearings considering striations 
running transverse to the direction of motion* Sxabsequently, 
Christensen and Tender [ 35 j derived a generalised form 
of stochastic Reynolds equation to include surface striations 
running in the direction of motion as well. A further refinement 
to reference [35] was made by assuming the flow fluxes to 
be represented by power series of a stochastic film 
thickness function- -Such a refinement "was proposed 
by Christensen, Shukla and Kumar . Sun [ 37 ]^ studied the 
effect of two dimensional surface roughness by taking 
into account the roughness spacing and autocorrelation 
function describing the roughness structure. 

The application of surface roughness theory depends 
on the characteristic length of the roughness structure. 

When this length is very much larger than the local film 
thickness the Reynolds equation is locally applicable ^ if 
;this length, is too small, Stokes equation is applicable [ 38 J. 
Sun and Chen [39 ]and Phan-thein [ 40 ] considered 
stokes roughness* Shukla [ 41 ] presented a deterministic 
approach to study the effects of roughness when the mean 
height of the surface asperities was of the same order of 
magnitude as that of the minimum film thickness* He 
introduced roughness interaction zone along the rough 
surface together with a purely hydrodynamic zone along 
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the smooth surface . This approach for roughness was 
later extended to the case where roughness interaction 
zones were considered along both the bearing surfaces 
and a purely hydrodynamic zone in between than [ 42 ], 

Recently, Pa tir and Cheng [ 43,44 ] presented an 
average flow model to determine the effects of three 
dimensional roughness by deriving the average Reynolds 
equation based on pressure and shear flow factors obtained 
through flow simulation. Brnploying this model and using 
asperity contact pressure obtained through Greenwood 
and Tipp' s model £ 45 Majumdar and Hamrock £ 46 ] 
studied the effects of roughness on the el as tohydro dynamic 
line contact by calculating the pressure for the entire 
contact region. Rec^tly, Prakash and Czichos £ 47 ] 
used Patir and Cheng' s model for partial el as tohydro dynamic 
lubrication of rough rollers. It has been pointed out 
that when there is no interaction of surface asperities, 
it is sufficient to analyse the inlet region,* however, 
in the mixed lubrication regime, it is necessary to analyse 
the entire contact region. Tender [ 48 ] calculated the 
average fluid flow and shear stress numerically in the 
hydrodynamic lubrication of surfaces consisting of adjacent 
periodic and symmetric unit roughness textures . 

Some experimental studies were also reported. Woo and Thomas[4 9]i 
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presented a review of experimental work relating to the 
area of real contact# the nxmber of contact spots# the 
spatial distribution of contact spots and the distribution 
of their sizes and the relation of all these to roughness 
and normal load. An interferometric study of the elastohy- 
drodynamic lubricated contact of rough surfaces was reported 
by Jackson and Cameron [ 50 J and scxne of Christensen's 
conclusions on the stochastic model for roughness [ 51 ] were 
also confirmed by it. A technique for measuring the heights 
of surface asperities in three dimensions has been developed 
by Tsukada and Sasajima [ 52 ]. It was confirmed that 
three dimensional characteristics of surface asperities 
on ground# lapped and worn surfaces could be studied by 
observing three-dimensional representation and contour 
maps obtained with the measuring system suggested by than. 

In the present wor]^ we employ Christensen' s 
stochastic model for roughness to study hydrodynamic and 
elas to hydrodynamic lubrication of certain conventional 
bearings using power law lubricant. 

1-4 elastohyderodynamic lubrication 

As mentioned earlier# under heavily loaded conditions# 
the bearing surfaces get deformed and the pressure generated 
in the film is modified due to elastic deformation and 
viscosity variation of the operating lubricant [ 53 ] . 
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The primary consideration in the el as tohydro dynamic (ehd) 
lubrication is the determination of minimum film thiclcaess. 
Grubin [ 54 ] was the first to study the ehd film thickness 
and obtained a formula for minimum film thickness. Dowson 
and Higginson [! 55 ] modified this formula after solving 
relevant equations numerically. However, it was pointed 
out that for all practical purposes these t^ formulae 
did not differ too much [ 56 ]. Laton on, Christensen [ 57 J 
modified the film thickness- formula for li^rication of 
rollers* Using Crook's approximation £ 58,59 ] and Grubin' s 
theory, many research workers have studied the ehd 
film thickness., Chow and Cheng [ 60 ] studied the effect 
of roughness with idealised asperities in the ehd contact 
of rollers. Dyson and Wilson [ 61 j investigated the 
problem of ehd film thickness of rollers, theoretically 
and experimentally using power law lubricants. Recently, 
Sinha et.al [ 62 ] analysed ehd lubrication of rough 
rolling surfaces with/without elastic deformation using 
power law fluids. 

In this work, we employ Grubin' s theory and Crook' s 
approximation to study the el as to hydrodynamic lubrication 
of smooth as well as rough rollers. 
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l-'o NON-NEWTONIiysI LUBRICATION 

As mentioned earlier# presence of high molecular 
weight additives (polymers) and dirt particles collected 
during manufacturing and lubricating processes make the 
operating lubricant in the bearing clearance deviate 
from the Newtonian behaviour. Such a lubricant is to 
be characterized by a non-Newtonian model £ 63-65 ]• 

Hirst and Moore [ 66 ] pointed out that many base oils were 
non-Newtonian tinder severe operating conditions. Hydrodynamic 
lubrication of gears and rollers involve rapid changes in 
the shear rate of lubricant and the relation between 
shear-stress and shear rate is no longer linear. In 
analysing the ehd film thickness, Johnson and Caneron £ 67 ] have 
pointed out that the fluids in high pressure ehd contact 
are essentially non-Newtonian in behaviour. During the 
last two decades several models such as power law, Binghan 
plastic Herschel-Bulkey, Reiner-Rivlin, pseudoplastic, 
viscoelastic etc. [ 68-74 ] have been proposed to characterize 
non-Newtonian fluids. 

In recent years, the power alw model has gained 
attention due to its capacity to characterize many types 
of lubricants such as polymer solutions and silicone fluids 
[ 61 ], Ng and Saibel [ 75 } -employed the modified 
model of pseudoplastic power law lubricant in the lubrication 
of inclined slider bearing and showed that the load capacity was 
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less than that for Newtonian fluids. Tanner £ 76 ] studied 
the non-Newtonian lubrication of short journal bearing using 
power law lubricants and showed that a reduction in the coeff- 
^ient of friction might be be obtained with typical power 
law fluids provided bearing eccentricity was not a -•limiting 
factor. Several attempts have been made to study lubrication 
problemsfor various bearing systems such as squeeze films [77 
externally pressurized bearings [ 78 ], conical bearings [ 79]^ 
journal bearings [SO], roller bearings [ 64 ] etc. Recently^ 

Kayer et«al [ 81/82 ] studied the behaviour of externally 
pressurized conical bearings with a power law fluid and 
pointed out that -the load capacity and temperature rise of 
the fluid increased wi’th the flow behaviour index of the 
fluid. In this work, the non-Newtonian behaviour of "the 
lubricant is characterized by power law model and a generalised 
Reynolds equation is derived and applied. 

\ -6 VISCOSITY/CONSISTENCY VARIATION I 

It has already been mentioned that/ for estimating I 

■the bearing characteris-tics in the actual performance of I 

the bearing system, generalisations to the Reynolds i 

equation to include effects such as viscosi-ty/consistency 
varia-tion, thermal effects etc. are to be made. Lope [ 83 ] modifiec’ 
the Reynolds equation by including viscosity and density 
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variation along the fluid film thickness. The viscosity 
variation along the film thickness has been considered 
by Zinkiewicz 84/85 ] . Dowson 86 ] generalized 
the Reynolds equation by considering the variation of 
fluid properties across and along the fluid film. Crook £87] 
assumed that the viscosity varied with pressure and 
temperature exponentially in the analysis of rollers. 

Thermal effects in lubrication have been recognized 
as an important factor affecting the bearing performances. 

A complete solution to the Reynolds equation was first 
obtained by Zienkiewicz by considering temperature variation 
across the film thickness with temperature dependent viscos- 
ity [ 84 J.McCallion et» al [88] studied the thermal 
effects by obtaining a complete solution to momentum 
and energy equations. Thermal effects are usually considered 
by assuming the lubricant to be Newtonian. However/ as 
mentioned earlier , due to a variety of causes the 

lubricant behaviour is essentially non-Newtonian. Hence, 
there is a need to study the thermal effects in the 
lubrication situation, using non-Newtonian lubricants. 

Many workers studied lubrication problems by coupling 
the Reynolds and energy equations. Rigorous consideration 
of temperature variation restricts the application of the 
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solution to a particular set of operational conditions 
i.e^/ to a particular set of viscosity tenperature 
relationship and to a specific inlet or initial viscrosity. 
Rigorous variable viscosity analysis is both expensive 
and limited in application [89]. Further, one of the 
major difficulties in calculating the tenperature of 
the bearing surface is that the number and size of the 
areas of contact between, bearing surfaces are not -usually 
known. Another difficulty is the experimental one of 
measuring the temperature at the contact zone in order 
to verify the proposed theoretical analysis [90]. 

The complexity of the complete solution to Reynolds 
equation mainly arises due to the inclusion of energy 
equation in the calculation procedure. Qvale and 
Wiltshire [9l] eliminated the energy equation from calculation 
procedure, by considering a chosen viscosity profile as 
an independent variable. This simplifies the calculation 
substantially and an extensive parametric study of 
lubrication problems can be attempted. A different 
approach to avoid the energy equation was presented by 
Tipei and his coworkers [92-94] who assumed a relationship 
between viscosity and film thicloness under thermal equilibrium. 
King et. al [95] remarked that a uniform temperature across 
the fluid film appeared to be reasonable in view of very small 
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film thiclmess C05~l4m) and very high linear velocity ^/s ) * 
Most of the heat generated in the film would be carried 
away by convection while thermal conduction across the 
boundaries would play a minor role» Following the 
approach of Tipei and his coworkers, in this work, thermal 
effects are considered by replacing consistency-temperature 
relationship by consistency film thickness relationship 
by a model akin to one proposed by them. 

in IN THIS WORK 

In this work, an attempt is made to study the non- 
Newtonian behaviour of the lubricant with the power law 
model* While using this model /many research workers 
derived generalized Reynolds equation without mentioning 
the domain for which it exists. Further, the non-dimensional 
scheme adopted by them are subject to controversy as the 
variables representing the viscometric indices of power 
law lubricant are used by them for non-dimen sionalisation. 
Obviating these controversies, Sinha et.al £ 64 ] derived 
a generalized Reynolds equation by assuming symmetry of 
flow of power law lubricant for lubrication of a cylinder 
on a plane. Prasad [96] considered rheological anomalies 
of the operating lubricant in the proximity of the bounding 
surface by considering the consistency of the lubricant 
film in the peripheral region different from that of 
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the film in the central region of the lubricant flow. 

The results were obtained by defining the load ratio 
as the ratio between loads with and without consideration 
of cavitation. Consideration of this load ratio 
parameter . . may lead to erraneous conclusions . 


Motivated by the need to define the flow domain 
for the application of Reynolds equation for power law 
fluids and to study the effect of consistency variation 
across the film thickness in the proper perspective, 
in this work we enploy the approach presented in the 
references [ 64,96 ], Consistency variation across 
the lubricant film is accounted through a general model 
for consistency variation which includes the variation 
of consistency due to thermal effects under the assumption 
of thermal equilibrium between the lubricant film and 
bounding bearing surfaces. This model [ 92-94 ] 
replaces the viscosity-tanperature relationship by the 
viscosity-film thickness relationship. Such a replacement 
is possible in tribological conditions as it has been 
experimentally verified that the highest temperatures 
occur in zones where the film thickness is least £ 92 ]■. 
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The model proposed in this work enables us 
to study the effects of consistency variation across 
as well as along the film thickness. It has been onployed 
to study problems of combined rolling and normal motion 
and el as tohydro dynamic Ixabrication of rollers. The model 
finds its application for smooth as well as rough bearing 
surfaces in the subsequent Chapters. Though the 
Chapters progress with a unified thene/ all of them 
are pres^ted almost as independent entities except for 
direct reference to equations derived in the earlier 
Chapters . 
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CHAPTER 2 


GENERALIZED REYNOLDS EQUATION FOR POWER LAW 
LUBRICANTS WITH CONSISTENCY VARIATION: 

APPLICATION TO ROLLER BEARINGS UNDER LIGHTLY . 

AND HEAVILY LOADED CONDITIONS 

2-1 INIRODUCTION 

Presence of additives in lubricants and consideration 

of high temperatures and pressures during bearing running-in 

give rise to variation of viscosity in the lubricating film 

across as well as along the film thictaiess [,1-3] • Keeping 

this in view,^ several atterrpts have been made in the past two 

decades to generalize the Reynolds equation considering 

viscosity variation* Dowson [ 4 ] generalized this equation 

by considering the variation of fluid characteristics (viscosity 

and density) across as well as along the film thickness* Ezzat 

and Rhode £ 5 ] accounted viscosity variation through a viscosity 

profile considering thermal effects. There are experimental 

evidences to support the need to consider viscosity variation 

in the fluid film under certain conditions. For example, 

Askwith et. al [ 6 ] reported that the presence of additives in 

the base oil led to a formation of rigid plastic layer in the 

proximity of the bounding surface in elastohydrodynamic lubrication 

The formation of high viscous layer very near to bearing surfaces 

has also been reported by Cameron and Gohar £ 1 ]•. Dyson pointed 

out that the thickness of the eihanced viscous layer might be of 
~8 —5 

the order, of 10 -10 cm* and ccxisideration of zones of enhanced 



35 


viscosity would have important effects in the lubrication of 
practical machinery [ 7 ], 

It is noted that the study of viscosity variation is 
generally carried out for Newtonian lubricants J scant 
attention has been given to this study for the case of non- 
Newtonian liibricants.Etasad [s] has enployed power law model 
for the liibricant with additives and considered consistoacy 
variation across the film thicloiess by assuming that the 
consistency of the fluid film in the peripheral layer 
adjacent to the bearing surface is different frcan that of 
the central layer. Little attention has been givei to study 
the effects of consistency variation which might arise due 
to adsorption of additives and variation of tanperature* 

Although a bearing performance is adequately predictable isothemaHy 

\ 

in many cases, the experiments showed consistent descrepancies , 
however small, when correlated with the existing theories* 

These theories should be modified to account for thermal 
effects and flow conditions [ 9 ]. King et. al pointed 
out that a uniform temperature across the film appeared to I 

be reasonable when small thicknesses and high velocities were 1 
involved. Most of the heat generated would be carried away 
by the lubricant while thermal conduction across the boundaries 
would play a minor role £ 10 ] . Keeping this in view, in this 
work, thermal equilibrium is assumed. Thermal effects, then, 
are accounted by assuming a relationship between consistency 
and film thickness akin to one proposed by Tipei £ 11 j.&uch 


36 


a relationship enables us to account for consistency 
variation along the film thickness in the undergoing 
analysis* 

Thus, consideration of additives and thermal effects 
leads to consistency variation across as well as along the 
film thickness- Keeping this in view^ in this chapter, a 
generalized form of Reynolds equation is derived. S\ifa.sequently, 
the effects of consistency variation on the iiabricatlon 
characteristics of roller bearings are studied- These effects 
are studied in the case of rigid rollers under lightly loaded 
conditions in S ec 2 . 3 ,* in Sec / the analysis has been 
extended to the regime of elas to hydrodynamic lubrication 
to study these effects on the minimum film thickness. 


2.2 GENERALIZED FORM OF REYNOLDS EQUATION 

Consider the symmetrical flow of a power law lubricant 
between two identical rollers, each having a rolling velocity 
U (Fig* 2.1). The thickness, h, of the fluid film in the 
lubricated contact is small compared to the radius r of 
the rollers. With the usual assumptions of lubrication theory, 
the basic eqns. governing the flow of a power law lubricant 

(in one dimensional foamn ) are given by [ 12 j 


^ = 
dx 


d . 

h 


i 6n 


1 1 
' dy ‘ 5y ^ 


^4., ^ ay., = 

ax dy 


(2 1 ) 


o 


( 2 . 2 ) 
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where u and v are velocities along the x and y directions, 
p = p(x) is hydrodynamic pressure and m = m(x,y) and n are 
viscometric parameters called consistency and flow behaviour 
indices, respectively* The lubricant behaviour is called 
pseudoplastic for n < 1, dilatant for n > 1 and Newtonian 
for n = 1. 

As the system is symmetrical about x axis, it is 
sufficient to consider the fluid region y ^ 0 to determine 
velocities. Taking note of the symmetry, the following 
boundary conditions may be prescribed for u ! 


au 

9y 


= O 


u = U 


at y = O 
at y = h/2 


(2.3) 


To determine the velocity u from eqn. (2*1), appropriate 
sign is to be attached to the velocity gradioit • To 

facilitate the determination of its sign, the velocity and 
pressure profiles are shown qualitatively in Fig. 2*2. 

It may be noted that the pressure becomes anbient at a point 
sufficiently away from the contact zone and reaches its 
maximxjm at a point x = -x* • In other words, 

^ =0 at X = -X* ( 2 * 4 ) 

Thus, in the region - oo< x<-x*, ^ ^ O and 3 ^ i 0 

^ dr) 

and in the region -x* < x < x* f *1 < 0 and 3 ^ < O, where 

— — G dy dx — 

x_ is the point at which the film starts cavitating. Thus, 


r 





ig. 2, 2 Lubsication of two rollers with 
law lubricant . 
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the velocity distribution for the region y > O can be 
obtained as 

u = -oo< x< -X* (2-5) 

and 

u = U+{- f "^y -x*<_x £ x^ (2.6) 

where the pressure p is denoted by p^^ and P 2 in the regions 
defined in eqpis. (2-5) and (2-6) respectively. 


Now, integrating the equation of continuity (2.2) 
with the conditions 

V = O at y = O 

V * Vj ^^2 at y = h/2 


(2.7) 


we have. 


h/2 


-V, 


h/2 


-9 r ^ U 
= gj- J u dy - ^ 


2 dx 


( 2 . 8 ) 


Note that V^y 2 resultant of the normal velocity V of 

the rollers and the normal velocity due to wedge action 
in the film segment, 


V - V 4. ^ 

^h/2 ~ ^ 2 dx 


(2-9) 


Positive values of V denote the normal separation of rollers 
while the negative values signify the normal approach implying 
squeezing action. 
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Substituting the expressions for u from eqns. (2,5) 
and (2*6) in eqn* (2.8), we obtain the basic eqns* determining 
pressure ' 

dx- -co<x<-x* (2.10) 

ix- I -X* < X < x^ (2.11) 

where 

h/2 , , 

(f) = f y(^) ' ^ dy {2.»12) 

o ^ 

In order to evaluate the integral on the RHS of eqn. (2.12) 
we have to determine the function m = m(x,y). 

It has been mentioned earlier that the consistency 
of the lubricant layer adjacent to the bearing surface 
may be different from that of the central layer and it may 
vary with temperature as well* The viscosity of all 
liquids, especially hydrocarbon lubricating oils decreases 
rapidly with increase in temperature* This variation in 
viscosity with tenperature is of practical importance in 
the lubrication of many mechanical devices such as gears, 
cams etc. where the lubricants are to function over a 
wide range of temper atures[ 13 ]• 

There is no fundamental mathematical relationship to 
predict accurately the variation in viscosity with temperature. 
The available viscosity- temperature relations are piirely 



42 


enpirical and the actual calculations require experimental 
data [ 14 3* In this study, it is assumed that thermal 
equilibrium exists and the consistency varies according 
to a given law. To apply this law to real lubrication 
situations, the temperature at each point should be known- 
This requires a complete thermal calculation. However, a 
viscosity- temperature relationship can be replaced by a 
viscosity-film thickness relation.- Such a replacement is 
possible in tribological conditions as it has been 
ejqjerimentally verified that the highest tsnperatures occur 
in zones where the film thickness is least £ 11 ]• In 
view of the viscosity-film thickness relation suggested 
in reference 11 ] , we assume that the consistency m of 
the power law lubricant varies according to the following 
1 aw ♦ 

— /h\ q 

" ’^^h^ (2-13) 

where in is the consistency at the inlet film thickness h^ 
where q, the thermal factor is measured, 

■q lies usually between O and 1. The value of q 
depends on the lubricant, the velocity of the bodies in 
contact, the cooling conditions and the flow of lubricant 
between bearing surfaces [ 11 Thus, q can be determined 

by completing the thermal calculations^ However, in the entire 
study of this work, various values are prescribed for q 
while examining the effects of consistency variation without 
going into details of thermal calculations. 
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It -will be recalled that the rheological behavioxir 
of the Ixibricant in the proximity of the solid surface 
results in an enhanced viscosity in the peripheral region. 
Particularly this phenomenon is evident when the additives 
in the lubricant are surface active and when the operating 
conditions of the lubricated contacts are in the elastohydro— 
dynamic regime £ 7 ]. This results in consistency variatican 
across the film thickness. Keeping this in view/- Prasad [ 8 ] 
considered the consistency variation by assigning that the 
consistency varies as a step function across the film 
thickness. Taking these two aspects# viz. consistency 
variation along as well as across the film thickness into 
consideration, we define, now, consistency as. 


m = 


= ">1 
1 


= Km, 1 

1 hi 


o < y < Y" ^ 




(2.14) 


where K is the consistency ratio and a is the thickness 
of the peripheral layer. 

Now, using eqn, (2,14) in eqn, (2,12) we obtain. 


(f) 


_ ^ 1 sl/n 

” 2n+i m^ 


/I'v (2n+l) /n 
^2^ 


j^q/n j^(2n+l-q)/n(f y 

1 o 


(2, 15) 


where 

(f^) = { l,.(l_2a/h) 


(2., 16) 
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The value of (f) when sxibstdLtuted in eqns* (2*10) and 
(2^11) yields the generalized Reynolds e< 5 i* which accounts 
for consistency variation I 

d r n _ ,1) (2n+l)/n , q/n , (2n+l-q)/n / . ) /I '^^l%l/n-| _ U ^ ^ ^ 

dxL2n+l ^1 ^ '‘^o^ ^itij^dx ^ ^ ~ 2 dx- ^ . 

— oo ^ fY> 0 (2.17): 

d (1) ^2n+l)/n , q/n , (2n+l-q)/n,. ) {_ 1_ ^^2a/n-, ' 

dxl- 2n»-l4^ ^1 ^ dx ^ -1 ! 

^ S (2.18),; 

i 

S 

2.3 LUBRICATION OP LIGHTLY LOADED ROLLERS UNDER ROLLING 

In this section, we study the effects of consistency 
variation on the bearing characteristics in the lubrication 
of two lightly loaded identical rollers under pure rolling 
conditions, with a power law lubricant, considering cavitation. 

The configuration of the system is as shown in Fig* 2.1. 

The Reynolds equation in this case, is written from eqns. (2.17) 
and (2.18) (putting V =0) I 

d r n /1% (2n+l)/n >,<3/n , (2n+l-q)/n wl ‘^^I'vl/n -i U dh 

d^ I 2n+i^F ' h (f^) ^ ' J - 2 d^ 

.^^x<_ -3^ (2*19) 

d r n ,lN(2n+l)/n v,q/n , (2n+l-q)/n/^ w 1 ^^2\l/n-i U dh 

2n+l%^ ^1 ^ W > - 2 d^ 

-X < X <_ Xq 


( 2 . 20 ) 
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To solve for pressure from eqns. (2.19) and (2.20), 
we use the following boundary conditions ! 


Pj^ = 0 at X = — oo 

( 2 . 21 ) 

dpj^ dpj 

dx ■ = ax = 0 at X = -x 

(2,22) 

Pj^ = P 2 at X = -X* 

(2,23) 

dp 

P 2 “ a'x = ° at X = x^. 

( 2 , 24 ) 


Condition (2«2l) indicates that the pressure is ambient 
at, a distance sufficiently away from the contact zone in 
the inlet region where pressure build-up starts. This is 
dependent on the oil quantity supplied per unit time» When 
the oil flow is above a required level, the film starts 
sufficiently away frcxn the point of contact* This point is 
virtually at infinity [ 15 ]. Eqns- (2.22) and (2-23) 
prescribe the matching conditions for pressure gradients 
and pressures, respectively, at x = -x* , the point of 
maximum pressure -Condition (2.24) locates the point at which 
the film steirts cavitating. In the general lubrication practice, 
it is assumed that the pressure is positive in the lubricated 
region and terminates at a point with zero pressure gradient 
where the film, forms a discontinuous mixtirre of air, vapour 
and lubricant (in the cavitated region )£ 14 ]- At the point 
of cavitation, the continuity condition of zero pressure 
gradient is satisfied as has been shown by Floberg 15 3* 
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(2.24) represents usual cavitation boundary conditions 
applied extensively in bearing analysis. 


Now, integrating eqns. (2.19) and (2.20) with the 
help of condition (2.22), we obtain the integrated form 
of the Reynolds equation as follows Z 


dp^ 

dx 


/2n+l 

^'2n“ 


. n 

) ■ 2 


2n+l 


m. 


hi 


(h-h*)"" 

’(£ h^n+l-q 

o 


--00 < X < wx’7y> O 
(2.25) 


fa = .(?n+l)n n <1 jq 

dx 2n 1 h^ 


( jnj^2n+l-q 
o 


X X^ , y> o 
(2.26) 


where the film thickness h is determined by approximating the. 
rollers by a set of parabolic cylinders at the point of minimim 
film thickness and is given by 

h = h^ + xV(2R) (2.27) 

where R = r/2 is the radius of the equivalent roller. Using 
eqns. (2.24) and (2.27) in eqn- (2.26), it can be seen that. 


x^ = X* (2.28) 

Thus, the points of maximum pressure and cavitation are 
synmetrical about the point of minimum film thickness. 

Integrating eqns.. (2,25) and (2,26) using conditions 
(2,21), ( 2 . 23 ) and (2.24), we obtain the expressions determining 
pressure I: 

I ' 
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Pi (x) = 


Po (x) = 


/2n+l\n m ttH/I 

2n 


y> 0 


-OO < X< -X 


(2.29) 


^2n+l\n -2n+l 
( „ ') 2 m 

2n 


7* -- I ^ a.] 

1 -OO jn^2n+l-q (f )nj^2n+l-q 

O ;• O 


-x*'<x<xt Y>0 (2,30) 

In order to calculate the point of cavitation x*, we use condition 
( 2 . 24 ) in eqn. (2.30) which yields. 


-X 


J l2(n,K,^q)dx = J l2(n,K,q)dx 


(2,31) 


-00 


-X 


where 


ll(n,K, q) = 


(H-H*)^ 


(P )“h 2"+1-<I ' ” 7“"n'H2n+l-4 

o o 


(2.32) 

( 2 . 33 ) 


and (Fq)= 1-(1-k"^/^) I l-(l-2i/H) J 

X = x/i2Rh^ , X* = x*/f2RhJ, H = h/h^ = 1+X^, H *= h^/h^ = 1+X 


a = a/h^ 

The solution of eqn. (2.3l) gives the point of cavitation- 


( 2 . 34 ) 


The load component W^(n,K, q) per unit loigth (in the x 


direction) on the cylinder is given by. 


W_(n,K,q) =j' ^ dx 


-00 


( 2 . 35 ) 
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Using eqns^ (2*.25) and (2*26) in eqni* (2.35) we obtain/ 


W^(n,K,q) = (Y-J" rr,, "r v ■ 

X ' 2n 1 R ^ J .jn^2n+l-q 

-oo o 


dx 


- r x" dx] (2.36) 

J tf )“ >,2n+l-q J 


-X 


(fo)“ h 


The load component Wy(n/K/q) in the y direction is given by. 


X 


W (n,K/q) = f P ” J* ^ ^ 


dx 


(2.37) 


-oo -oo 

which on using eqns. (2.25) and (2.26)yields, 


W. 


,(n,K,q) = -(^5;!:^)" ’/X -'■fi-i'-*'." • • - 

^ *'1 •' (f ,n ^2n+l-q 

o 


dx 


-oo 


(2.38) 

To Study the effects of consistency variation across the 
film thichness on load, we define the quantities and 

as follows “* „ 


-X 


X 


W (n 


w. 


V / X^I, (n,K,q)dX- / X^I„(n/K,q)dX 
_ 0 O -X 


Kx ■ w^(n;r,q) „ 

^ -X 


(2.39) 


X 


X I^(n,K,q)dX- / X^^l 2 (n, K, q) dX 

-X* 


«»oo 


and 


W^(n,K/q) 
^KY- ^ ^Wy(n,l, q) 


^K 1 

-X* 

/■ XIj^(n,K,q)dX 


-oo 


(2.40) I 


“^1 

f ^ XI^(n,l,q)dX 
-oo 



49 


where the expressions l 2 _(n/K,q) and l 2 (n,K,q) are defined in eqn; 
(2*32). 


The case K = 1 indicates no variation of consistency across 
the fluid film J that is, the absence of peripheral layer in 
the lubricant film* It gives the consistency variation along 
the film thickness only* The points and denote the 

points of maximum pressure and cavitation, respectively, for the 
case k = 1* 

The effect of consistency variation along the film thickness 


on load can be studied by defining the quantities and 


as follows ! 

W^(n,K,q) 


W^(n,K,07 


* * 

—X X 

J X^I^(n,K,q)/Hj dX - J X^l 2 (n,K,q),/Hj dX 


-oo 


-X 

■ 

/ ° X^I^(n,K,0)dX- / ° X^I (n,K,0)dx 


(2*41) 


-X, 


-OD 


-X 


qo 


-X 


W ^ = 
qY 


W^(n,K,q) 

Wy(n,K,0) 




f X I^(n,K,q)/Hj dX 


(2.42) 


-X 


% 


f XIj^(n,K,0)dx 
-oo 

The case q = O corresponds to the isothermal consideration? and H^=h-j. 
Thus/it gives the consistency variation 'across the film thickness oni 
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The points and represent the. points of maximum 

pressure and cavitation, respectively# for the case q = O- 
The e^qaression for frictional force on the surface y = h/2 
is given by# 


X 


fo = /■ T 


-oo ly=h/2 


dx 


(2-43) 


where 


. I ftv I av 


(2 ,44) 


dy ‘ dy 
T, being the shear stress. 

Using the expressions for u from eqns* (2*5), (2,6) and (2.44) 
in eqn. (2.415), we obtain 


, dp. . dp 

^2 “ “I / I dx^ <lx + / ^ dx] 


2 dx 


(2*45) 


-oo -X*" 

With the help of eqn. (2.35) it can be seen that# 


f 2 (n,K,q)-- W^(n,K,q)/4' 


■ (2-46) 


Similarly, the frictional force fj^(n, K#q) on the surface 
y = h/2 is given by, 

f 2 (n,K,q)=-W^(n,K, q)/4 (2.47) 

Thus we have# 

f 2 ^(n#K,q) + f 2 (n,lC#q) + W^(n,K,q)/2 = O 


(2.48) 
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This is the condition for overall equilibrixm that must be 
applied to the lubrication of roller bearings in the absence 
of inertia and body forces- Siibstituting in eqn- (2-48) 
k = 1 and q = O we obtain the. result of reference [16] 
in the case of Newtonian lubricants^ 

2.4 ELASTOHYDRODYNAMIC (EHD) LUBRICATION 

In this section/ we study the effects of consistency 
variation on the minimum film thickness in the ehd lubrication 
of heavily loaded rollers. The analysis is restricted to 
the inlet region. For this purpose, we use Grubint theory [,17] 
and Crook's approximation [ 18 ] in the analysis. In Griobins 
theory tvro simplifying assumptions are made [ 16 ]!• (1) 

The deformed shape of the bearing surface is the same as that 
in the dry contact- (2) A high pressure is developed in the 
entry region to the main loading region, called Hertzian 
zone- The physical situation is as shown in Figv2-3 where 
2d is the width of the Hertzian region and the pressure is 
maximum just outside this zone, i.e., at x = -d. 

The eqn. to determine pressure Pg in the inlet zone of 
ehd regime is written from eqn- (2.19) with p^^ replaced 
t>y Pg as. 




rn^(h /h|)^ 


y 

A 


/ . 


X 



xs-d Hertzian zone 


Fi< 2.3 Eiastohydrodynamic lubrication betweGn two 

rollers considering consistencyvariafion/!,j..w;’a 
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dr n (2n+l)/n . (2n+l-q)/n ,r wl '^^E^l/nr U dh 

&'2n+l ¥ ^ ‘S7 rx'^ ' ] = 2 K 


-oo ^ X ^ -d ( 2 . 49 ) 

where £ 15 ] 

d = 2[2WgRg/( 1/E = (l-^)/E^, Rg = r^2 (2.50) 


y is the Poisson's ratio, E and E^ are the Yoimg' s moduli.., 
Rg is the radius of the equivalent roller, r^ is the radius 
of the lubricated roller and Wg is the load in the ehd case. 
To determine pressure in the inlet region from eqn. (2-49) 
we use the following boundary conditions I 


Pg = 0 at X = -00 (2.51) 

dp„ 

Pg = 00 , = 0 at X = — d ( 2 . 52 ) 

Integrating eqn, (2.49) with the second condition of (2,52), 
we obtain 



dp. 

dx 


E , 22"+1 ra^u" ■ - 


2n 


(h-h ) 
m 


n 


( 2 . 53 ) 


(f )^ i^2n+l-q 


where h = is the minimum film thickness at x = -d which 
is constant along the Hertzian zone. To write an expression 
for the film thickness in the heavily loaded lubricated contact, 
we first note that the deformed shape of the rollers outside 
the band of Hertzian contact is given by | 20 ] 

h2 = (d V ( 2R))[ i 3^d [ ( x^/d^-1) ^/^-In C I 1 + (x^/d^ -1 ) } ] 

(2.54) 
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which may be approximated as [ 18 ] 

h 3 = 2^/^ d^(3R)"^ + higher powers of e (2.55) 

where e is given by 

I x[ = dCl+e) (2.56) 

According to Grubin' s theory# in the Hertzian region the 
bearing surfaces are separated by a constant film thickness. 
Hence# the film thickness in the ^d contact of rollers is 
given by the deformed shape of the rollers outside this 
region with the addition of constant separation in the 
Hertzian region. Thus# we have 

h = 2^/^ d^CSR)"^ + h (2.57) 

m 

In the ehd lubrication the consistency of the lubricant 
is assumed to vary exponentially with pressure [20# 21]! 

mj^ = m^ e°''^ (2.58) 

where ct is the pressure coefficient of consistency index. m^^.This 
type of eqn. is beleived to give a fairly accurate viscosity- 
pressure relation in the case of high viscosity fluids under 
isothermal considerations [ 22 ]. The consideration of 
consistency variation with pressure under isothemal conditions i 
eqn. (2.58) is valid# for# in the definition of m in eqn. (2.14) 
m^^ may principally be a function of pressure £ 11 ]. 

Now# integrating eqn. (2»53) using eqns. (2 .54) -(2.58) 
we get# 
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i = 2^"+! n. U” . ax (2.59) 

° *^1 -oo (f ,n j^2n+l-q 


a 


(2.60) 


On simplifying, the formula for minimum film thickness in 
the case of ehd lubrication can be written as follows ! 
\)Vjn:rx;/j ^ q 

where 

, C, = (2n«)" 23/2 3 -V3 ^1/6 ^ 


(^)(3rH-l)/3 ^ (n,K.,q,0) 


lt/2 


m 

(2.61) 


I^,(n,k,q,e) = / (cosS)'®"-®'5-3V3!(3i„9)(6n+l)/3/(j, ) 

jiX 

(2.62). 

(P ) = l-Cl-K"^/’^) £ l-(l-Pcos^e ) (2.63) 

o 

h = h sec^e , P = ^ / d^(2e)^/^ = tan^e(2.64) 


n 


m 


Eqn- (2. 60) gives an expression for the minimum film thickness 
in the case of 1^^ 2a. The minimum thickness formula for the 
case of no consistency variation across the film thickness 
(k = l) is written from eqn. (2*60) as. 


h^(K-l) (3n+i)/3 


= C, am^ (^) ^ (n, q, %/2) (2 . 65 ) 


■E2 


where 


Ig 2 (n,q,e) = / (cose)^^’^"®'^"'^^/^' (sine ) ^ de 

° ( 2 . 66 ) 


Using eqns. (2.60) and (2.65) we write the ratio 
in the case of h^^ ^ 2a as. 


m 


h 


m (k=l) 


de 
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^(k=1) “ (2-67) 

m 

The ejq^ression for minimim thiclcness formula with no consistency 
variation along the film thickness is given by l^(q=0). Hence, 
we write the ratio h^/h^(q=0) for as 

h^/\i^q=0^ =[(1/H^)'^ Ig^(n,K,q,0)/Ig^(n,K,0,0)]^^^^”+^^ 

( 2 . 68 ) 

When h^l,2a, the inlet zone is divided into two regions 
-CO < X <_ -dj^, h > 2a and -d^ < ^ 1 ^ 1 shown 

in Fig.2»4 where h = 2a at x = -dj^» 

The governing pressure eqn* for the region -oo <_ x <_ -d^^ 

is rewritten from eqn. (2. 49) by taking p. for p as 

XE E 

d ir n /l^(2ntl)/n q/n (2n+l-q)/n/^ wl_ ^^lEvl/n-i U ^ 
dxL2n+i ^1 ^ '■^o ra^ dx ^ 2 dx 

(2.69) 

In the region -d^^ < x <_ -d, the consistency of the fluid 
film varies as icm^^Ch/hj^)^ in the entire region. The eqn. 
for pressure P 2 g in this region can be written as 

d r n /l'i(2n+l)/n r^g/n , (2n+l-q)/n r 1 ^^2E\l/ni U ^ 

aSjL 2S.l(2> ' hj- h j=2te 

(2.70) 

The boundary conditions for eqns. (2.69) and (2,70) can be 
written as 
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Fig.2.4 Elastohydrodynamic lubrication between two rollers, 
considering consistency variation , hyyi ^ 2a. 
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P 2 _g = P at X = -cx) 


(2.72) 


^lE ^2E ^ ~ **% 


( 2 . 73 ) 


dx = P2E = X = -d 


( 2 . 74 ) 


Integrating eqn. (2.69) we get 


^lE _/2n+l\n ^2n+l TT^i/h %q 
dx 2n 1 hj^ 


X [-. - - ^ ? . . ..... f 

^(2n+l)/n__|(j^__j^-l/nj ^j^(2n+l) /n 2 g)^ 2 n+l) /n } 


(2.75) 


where is the constant of integration to be determined.. 


Integrating eqn, (2.70) and using condition (2.74) 


|e ,(2n*l,n 32 n.l ,, ^ jn 


where h = at x = -d. 


Now, using condition (2.7l) in eqns. (2.75) and (2.76) we get h =-h , 

ap, “ 

(Xp . 

Recollecting that m, = m e , we use m, = m e xn the 


xr^ ET /* 

region -co< x ^ and = m^e in the region -d^^ ^ ^ d. 


and integrating eqns. (2.75) and (2.7 6) using conditions 


(2*72)-(2.74) we get 


= aC 3 _m^(R/h^) ^3^+1 (n,K,q,e ),/H? (2.77) 


= Kac,m (R/h ){3n+l)/3 I (n,q,ei)/H? (2.78) 

1 o ^ rn E 2 * ^ 1 
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where I^Cn, K, q,G -j^) and are given by eqns., (2. 61), 

(2.62) and (2,66) respectively/, and 

= cos P = 2a/]:^* (2.7 9) 

Adding eqns, (2.77) and (2.78) we have 

(rl)(3n+l)/3 = aC^m^£lg^(n/K/q/8;j_^+^^E2^'^'^^®l^^/^l 

(2.80) 

Eqn. (2.80) gives the mini. VAtn thickness formula for h 5 2a. 

nr^ 

We define thu ratios h. /h (k= 1) and h /h (q=0) in this 

m' m m ^ 

case as well^ 



h^CK=l) 



(h/K/ q/6 
1^2 ^n/qrTc/2'^ 


(n/q, Sj^) 


j3/(3n+l) 

(2.81) 


h 

m 

h (q=0) 
m ^ 



Igj^ (u/K/q/ ©2^) + -^^22 ® 

Igl (n/Kv 0/ Sj^) "^^^22 ^1^ 


j3 /(Sn+l) 


(2.82) 

We see that putting q = 0 in eqns. (2,67) and (2.81.) .we 
obtain the minimum thickness formulae for h^ 2a and 
h„ £ 2a as shown in reference [ 8 ]. 
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2.5 RESULTS AND DISCUSSION 


Eqn- (2.31) is nijmerically solved to determine the 
point of cavitation, X'^ ^ for rigid rollers ^for various values 
of the flow behaviour index n, the consistency ratio K, 
the peripheral layer thickness a and the thermal factor q. 

The load ratio parameters ^KY' ^qX ^qY thei, 

calculated by using eqns, (2.39)- (2. 42) , The case K>1 
signifies an increase in the consistency of the lubricant 
in the peripheral region, increasing thereby the effective 
consistency in the lubrication process relative to the case 
with no peripheral layer. The case K< 1 indicates a decrease 
in the effective consistency. The latter situation may arise 
with certain additives when chemical reaction takes place 
near the bounding solid surface. This will cause the additive 
particles to disappear from the bounding surface or the 
lubricating film [23]. Davenport [ 24 ] pointed out that 
higher viscosity might be due to the presence of additives 
or surfactants in the carrier oil. This condition, in our 
study, pertains to the case of K> 1. 


Increasing value of q signifies a decrease in the consistency. 
This may be a consequence of temperature rise as it is well 
known that temperature rise results in viscosity/consistency 
reduction. The case q = o represents the isothermal consideration, 
and provides maximum pressure geieration and load capacity. 



Throughout this work, we represent themal effects by q, 
and call it as thermal factor. Setting q = o in the 
present analysis, we obtain the analysis of Prasad £ 8 ]. 
However, his results cannot be compared to that obtained 
in this study because the load ratios and defined 

by him carry dual consideration of cavitation and no 
cavitation. These ratios lead to erraneous conclusions 
for, consideration and disregarding of cavitation will 
alter the load capacity significantly. In the present 
analysis, the ratios are calculated considering cavitation 
only. The value of n normally lies between 0.'^ and 2.5 £25]. 
For the purpose of comparison of pseudoplastic and dilatant 
behaviours, the values of n =.0«5 and n = 1. 5 have been 
taken as the representative values. 

Fig. 2. 5 gives the picture of the location of cavitation 
point with respect to specified K for various values of 
n. As K increases, it can be observed that X* moves 
towards the point of minimum film thickness; however, 
the shifting towards the converging gap is very small for 
increasing values of large K. Further, for any fixed 
K, the cavitation point for ^illatants lies closer to the 
point of minim'um film thickness as compared to pseudoplastics 

From Fig* 2.6, we observe that for K< 1, an increase 
•• ^ 

in a results inashift^uja/ of X frcxn the point of minimum 
film thickness for all n ^ for K > l,the trend of shifting 
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is reversed. 

The effect of q on the location of X * is observed 
in Fig. 2. 7. An increase in the value of q shifts the 
point of cavitation towards the region of diverging 
film segment. As pointed out earlier, increasing 
q signifies reduction in the consistency and pushes the 
cavitation point towards the diverging gap. It is 
observed from the Figure that the effect of q on X* 
is more pronounced in the case of dilatants compared 
to pseudoplastics for large values of q. 

Figs. 2.8 and 2.9 depict the effect of K on load 
ratio parameters and ^^om these Figures, we 

observe that the effect of K on load capacity for 
dilatants is much more pronounced compared to that for 
pseudoplastics . This effect gives greater load 
capacity for K > 1 and less load capacity for K< 1 
for dilatants compared to pseudoplastics. The effect 
of increasing a in the peripheral layered film thickness ( K>l) 
provides greater load capacity compared to the case with 
no peripheral layer. A reversed trend of load capacity 
is observed in the case of K < 1 (Figs. 2.10,2.11). 

From Figs. 2.12 and 2.13, we observe that as q increases 
the load capacity decreases for all n. The decrease is more 
for dilatants ccanpared to pseudo pi as tics for large 
values of q. 
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The effect of consistency ratio K and thermal factor 
q on the minimum film thickness in the ehd lubrication 
are studied in Figs* 2.14~2.17. To study the effect 
of eqns. (2.67) and (2.81) are numerically evaluated 
for hj^2a > 1 and h^/2a < 1 respectively for specified 
values of h^/2a for K = 1. From Fig. 2.14 we observe 
that the effect of K is to increase the minimum film 
thickness. For K > 1, the film thickness will be higher 
than that of the case with no peripheral layered ehd' 
minimum film thickness and for K < 1 the film thickness 
is smaller. 

The effect of K on minimim film thickness for various 
values of n is manifested in Fig. 2. IB*. For K< 1, the 
increase in the minimum film thickness is more for 
pseudo plastics compared to that dilatants. This behaviorur 
is in qualitative agreement with the findings of Prasad [s]. 
The trend is reversed in the case of K < 1 when h^(K=l)/2a 
is sufficiently large. It is to be noted that no definite 
pattern of the flow behaviour is observed when the 
peripheral layer is close to the half of minimum film 
thiclcness in the case of K < 1. 

The effect of q on minimum film thickness is studied 
with eqns. (2.68) and (2.8'2-). An increase in the value 
of q for K> 1 reduces the minimum film thickness (Fig. 2,16) . 
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However for K <1, there is slight increase in the minimum film 
thichness(F(^2,17 ) * One may infer that for K ^ 1 an increase in the 
thermal factor has a 'beneficial effect on the ehd minimum film 
thickness. The calculations showed that in the case of ehd lubri- 
cation, a reasonal prediction of film thickness is not always 
possible with high values of q. Hence the formulae for minimum 
film thickness will be reliable for low values of q, particularly 
for isothermal case. 

2.6 CONCLUSIONS 

In this Chapter, a generalized one dimensiohal Reynolds equat- 
ion is derived for power law lubricants with consistency varying 
across and along the film thickness. Consideration of consistency 
variation across the film thickness is due to the presence of 
peripheral layer and this variation along the film thickness 
is due to thermal effects. 

The effect of the high consistency peripheral layer is to 
shift the cavitation point towards the point of minimum film thick- 
ness and to increase the load capacity with respect to flow behav- 
iour index. The trend is reversed for the case with low consistency 
peripheral layer. The effect of thermal effects is to decrease 
the load capacity,* the decrease is more for dilatants than for 
ps eudo pi as ti cs . 

For heavily loaded rollers, the effect of high consistaicy 
peripheral layer is to increase the minimum film thickness,* the 
increase is more for pseudo pi as tics compared to dilatants. Th,e 
trend is reversed for low peripheral layer film thickness, provided 
the film thickness is sufficieitly high* , 

t 
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No set pattern of the flow behaviour is observed when 
the peripheral layer is close to half of ^d film* The 
effect of thermal factor is to decrease the minimum film 
thicl<ness in the case of high consistency peripheral 
layer. The traid is reversed for low consistency peripheral 
layer. 
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Fig. 2.6 Cavitation point X vs. peripheral layer thickness 
a for various values of flow behaviour index n. 
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Fig. 2.8 io'id capacity ratio vs. consistency ratio K for 
VO ous values of flov^ behaviour index n. 





Fig. 2.9. Load capacity ratio Wj^y consistency ratio K lor 
various values of flow behaviour index n. 
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Fig. 2.11 Load capacity ratios and Wgy vs. peripheral 
layer thickness 5 for various values of flow beha- 
viour index n, K>l. 
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NOMENCLATURE 


a 

a 

d 

% 

E 

h 

h* 

H, 


K 

n 

P 

P1/P2 


Pe 

Pm 


peripheral layer thickness 
non-dimensional quantity correspond to a 
half width of the Hertzian contact 
point on the inlet region where h = 2a 
reduced modulus of the two cylinders - 
Young^ s modulus of the cylinders 
frictional forces at y = -h/2 and y = h/2 
respectively 
film thickness 

minimum film thickness in the hydrodynamic 

lubrication 

inlet film thickness 

film thickness at the point where pressure 
gradient is zero 

minimum film thickness in the ehd Ixibrication 

non-dimensional quantities corresponding to h, h^^ 

and h respectively 

consistency ratio 

consistency indices 

flow behaviour index 

hydrodynamic pressure 

pressure defined in the regions with positive 
and negative pressure gradients 
pressure in the ehd lubrication 

pressure at x = -d^ (eqn- 2,73) 
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PlE' ^2E 


r 


^E 

R 


R, 


U* V 

u 


V 


V, 


h/2 


\ y ' \y 


ehd pressures in the regions defined 
in eqns. (2- 69 ) and (2. 70 ) respectively 
radius of the roller 

radius of the roller in the ehd lubrication 
radius of the equivalent roller 
radius of the equivalent roller in the 
ehd lubrication 

velocities along the coordinate axes 
rolling velocity 
normal velocity 

resultant normal velocity at y = h/2 
X and y components of load 
non-dimensional load ratio parameters 


x,y 


-X 


x^ or x’^ 

-X*, X * 


* * 


T 


coordinate axes 

point of maximum pressure 

point of cavitation 

non-dimensional quantities corresponding 
to X, -X* and x* respectively 
non-dimensional points of cavitation for 
the cases K = 1 and q = 0 respectively 
shear stress 


V 


Poisson ratio 



CHA^^ER^ -3 

EFFECTS OF CONSISTENCY VARIATION OP 

POWER LAW LUBRICANTS IN SQUEEZE PII14S 

3.1 INTRODUCTION 

The study of squeeze films has been one of the subjects of 
interest in lubrication applications and rheological studies. 
However, such a study is not of great practical value, but it 
is representative of time dependent lubrication [ 1 
Theoretical work on Squeeze films using a Newtonian fluid 
on various geometric configurations was reported by 
Archibald [ 2 ] . Several investigators have analys-ed squeeze 
film lubrication using non-Newtonian fluids These 

studies/ in general, do. not take into consideration variation 
of viscosity/ consistency of the operating lubricant- Such a 
consideration is warranted by an increase in the effective 
viscosity, particularly in the lubricant layer very near 
to bearing surfaces. Needs [ 6 ] observed a predominant 
enhancement of viscosity dinring squeezing- To detect the 
influence of bounding surfaces on viscosity of the squeezed 
films, he conducted a series of experiments using two 
optically plane parallel circular plates approaching each 
other with a normal velocity, and measured film thickness 
down to 0.635E-03 mm. He noticed an increasing discrepancy 
betweei the measured and theoretical (as measured from 
classical Newtonian theory) intervals of time (response time). 
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for film thicknesses less than 0.127E-02 mm. The increase 
in the effective viscosity in such a thin film was 
attributed to the influence of metal surfaces on the 
fluid layer in their con vicinity, causing the fluid to 
behave more rigid. Derjaguin [ 7 ] observed an increase 
in viscosity in the liquids where polar surface- active 
substances were present at a distance upto O.lum from 
the solid boxindary. He pointed out that in liquid 
polymers and in polymer solutions, there is generally 
an increase in viscosity at a distance upto 7 or 8 M-m 
from solid boundary. Hayward and Isdale [ 8 ] attributed 
the abnormality in rheology to dirt and ranarked 
that in pure liquids this would not happai. In the 
experiraantal findings of Askwith et.al. 9 ] it has been 
pointed out that the organic liquid in contact with the 
bearing surface formed a high viscous layer adjacent to 
it. Yousif et.al [ 10 ] observed that the molecules of 
the liquid in intimate contact with or to the adsorbed 
layer of the bounding surface must behave differently 
from the adjacent molecules in the bulk lubricant. However, 
the entire nature of the surface of the influence and 
the distance to which it penetrates has not been fully 
understood. 

It is evident,on the whole, that there exists a strong 
case for the consideration of consistency variation in 
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squeeze films. The effects of consistency variation 
have been considered in Newtonian fluids [11-14]. Scant 
attention has been given to such a variation using 
non— Newtonian fluids [l5]. With particular reference 
to squeeze films, Shukla et.al [ 15] analysed the non- 
Newtonian behaviour of lubricants through power law model 
by considering consistency variation across the film 
thicl<ness* In this chapter, we consider the consistency 
variation across as well as along the film thickness by 
adopting the model we proposed for such a variation in 
Chapter 2. In particular, we study squeeze films on rigid 
solids for various geometries such as parallel plates, 
circular plates, roller and journal bearings etc. 

and the effect of consistency variation on load and response 
time is analysed in each case«- 


3.2 PARALLEL PLATES 

In this section, we consider the flow between two 
parallel plates of length 2'd, approaching each other 
normally with a velocity V due to a symmetrically placed 
load (Pig. 3.1). The plates are rs'eparated by 
a film thickiess 2h, With the usual assimptions of 
lubrication theory,' the governing eqns. of motion for a 
power law fluid in the case of squeezing is obtained 
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Fig. 3J Squeeze film between parallel plates. 
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with reference to Pig. 3.1 and putting U = 0 in eqn. (2.18). 

where 

(f^) = 1-(1-K~^/^) il-d-a/h)^^^"^^^/^ } (3.2) 

and h^ is the initial film thickness measured at x = -d 
just before squeezing commences. 

Pressure attains its maximum at x = 0, i.e., 

at X = O.Using this condition in the integration of eqn. (3.1) 
we get# 


t = -h '"n- '3.3) 

o 1 

Integrating again eqn. (3.3) using condition p = 0, at 
X = d# we obtain the expression for pressure p. Denoting 

it by Pj^ we have 

/ g 


m 


%q 


1. (?.n+l V ,n (1 ,g (l,2n+l-q (^^+l_^n+l) 


n+l n 


(fo) 


(3.4) 


The load capacity ^ yier unit width is given by 

d ■■ _ 

W„ = 2 / Pp. (x)dx (3.5) 

v/hich on using eqn. (3.4) becomes 

„ , - !”ll (?ntl V )n 311+2 (^,q (^,2n+l-q 

K^q n+2 n (f^) h^ h 

The squeezing time ^ from an initial film thickness 
2hQ^ to a subsequent film thickness 2 h 2 # say# is obtained 
by putting -V = in eqn. ( 3 . 6 ) and integrating# we have 
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= Sn+l (n+ 2 )/n ( jl/n / < 5 h 

n Cn+P^^--' ^ ^ 


'K,q “1 h 2 (f')i;^2n+l-.q)/n 

( 3 . 7 ) 

Considering consistency variation along the film thickness 
only (which can be obtained by taking K=l),eqns. ( 3 . 6 ) and 

(3.7) beccrae- 

2m- - _ _ _ - ^ 

( 3 . 8 ) 


1 / q n +2 n h^^ h 


and 


2 n+l ,(n+ 2 )/n _^1 a/n (l_)q/n 


n +2 W, 


'1, q n+l-q q 

X [ { 1 . ) - ^1_) (n+l-q) /n ^ (3^9) 


To determine the effect of thermal factor q on 
load and response time, we define the following 
quantities Z 


where 


11 

%Sl = h '^ 

^K,o 

(3.10) 

t = 

q 

^S/g= 

^,0 ^K,o 

(3 . 11) 

^K,q 


(3.12) 

(f^) 

0 

= l-d-K'V’^) { l>(l>i/H) 

(3.13) 

H = 

h/h 3 _, a = a/hq, = ^ 2/^1 

(3.14) 
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To study the effect of consistaicy ratio K on load 
and time of squeezing, we assume that the peripheral 
layer thickness varies with film thickness only, i. e. , 
a = 5 h where 6 < 1 is a constant, and define quantities 


K 


and tj^ from eqns. (3. 6) -(3 -9) as follows 


W, 


and 


W = 
K 




K 


W 




n 




V 




tq = -- 


'i/q 




where Fg-x= 1- (l-K'^H) £ (2rH-l)/n j 


(3.15) 


(3.16) 


(3.17) 


Eqn , ( 3 . 10 ) gives the variation of load ratio with q. 

Since W is a function of t and q and at any instant of 

q 

time during squeezing H< l,thus,W^ decreases as q 
increases* This is in conformity with the established 
result that the thermal effects tend to decrease the 
pressure generation and hence the load capacity. Eqn. (3*11) 
gives the variation of t with q for various values of n 
(Pig. 3 . 2 ) . As q increases the response time decreases for 
any fixed but arbitrary n, possibly due to the thinning 
of lubricant with thermal effects. These effects significan- 
tly alter response time for pseudoplastics compared to 
dilatants. 

To study the effects of K and a on and tj^,eqns, (3.15) 
and ( 3 . 16) are evaluated. The effect of K on is studied 
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in Fig,3*3whGre the peripheral layer thickness a varies in 
constant proportion to half of film thickness 2 h. The 
effect of an increase in consistency in peripheral region 
is to increase load capacity for all n. Further, load 
capacity is more altered for pseudoplastics compared 
to dilatants. From Fig, 3, 4 we note that the change in K 
alters the response time more for pseudo pi as tics than for 
dilatants. 

The effect of K on load and response time is further 
eloborated in Figs. 3.5 and . 3*6 respectively for specified 
value of h. For fixed h, increase in 6 * means an increase 
in peripheral thickness and load capacity and response time 
increase or decrease according as K ^ 1. Further, for 

specified film thickness and K, load capacity and response 
time are altered more for pseudo pi as tics compared to dilatants 
due to variations in peripheral layer thickness. 

3 .3 PARALLEL CIRCULAR PLATES 

We consider squeezing between two circular parallel 
plates each of radius R separated by a film thickness 2h 
(Fig. 3 -7). The plates approach each other with a normal 
velocity V symmetrically. The Reynolds ..eqnati.on for the flow 
of power law lubricant in the radial direction r can be 


obtained as 
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1 d r " (f ) r(- 1- |f)V"]=v 

r dr ^ 2n4-l 1 o 


(3.18) 


where (f^) is defined in eqn. (3.2). 
Using the boundary conditions 


^ = 
dr 


0 at r = O 
and p = O at r = R 

we can obtain an expression for pressure by integrating 
twice, eqn. (3.18). Denoting it by p^^^ ^ we have 

^ "^1 (2n+l V ^n (l_)q A) 2 n+l-q(Rn+l__^n+l) 

%,q n+1 ^ 2n Cf^) h^ h 


(3.19) 


The load capacity is given by 


p.K 
R 

W„ = J’27irpdr 

Q Q 


(3.21) 


which on using eqn. (3.20) yields. 


_ ’^■^l /2n+l V ^■^n j.n+3 (1_)<3 - 

^K.q ' n+3’‘ 2n If^) ^2n+l-q 


( 3 . 22 ) 


The time of squeezing g "to reduce the initial film 
thickness 2h^ to a subsequent film thickness 2 h 2 is given 

by 




2n-H , ’""l V )V" Rln+sj/naj-l/n 

- 2n '"n+S t>l 


X 


^1 . 
/ 1 


(f )h^2ntl-:q)/n 


dh 


(3*. 23) 


X 
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As in the case of parallel plates, the quantities to 
study the effects of q and K on load and time of squeezing 
can be defined as 


W, 


and 


W = 

q 




w 


= H- 


K,0 


t = 


'K,0 




K,o 

1 jn 


K Wi 

l#q F -jf 
o 

t. =. = -A- 

^ h,q ^6* 

where ^ and F^* are defined in eqns- (3.12) and 
(3.17) respectively. 


(3.24) 


(3.25) 


(3.26) 


It is seen that the expressions for bearing characteristics 
obtained for parallel circular plates (eqns. (3.24)-(3«27) ) 
are the same as those obtained for parallel plates ssqns. 

( 3 . 10 ), ( 3 . 11 ), ( 3 . 15 ) and ( 3 . 16)). Hence,, similar inter- 
pretations can be made in this case as well. 

3 .4 ROLLER BEARINGS 

The problem considered is that of the flow of a power 
law lubricant between two identical rollers each of radius 
r which approach each other with a velocity V (Fig. 3.8 )« 

The rollers have a projected length of fluid film 2d and 
are separated by a film thickness 2h. The film thickness 
is small compared with the radius of the rollers. With the 
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usual assumption of lubrication theory and. symmetry of 

bearing configuration, the Reynolds equation can be written 

in the case of squeezing for rollers using eqn. (2.18) and Rig * 3 , 8 !; 


d r n_ q/n ^( 2 n+l-q)/n 
dx*- 2 n+l '^l ^ 



(3.23) 


where 2 hj_ is the initial film thickness measured at x = d 
and h is given by 

h = h^ + xV(2R) (3.29) 

2h^ being the minimum film thickness and R = r/2 the radius 
of the equivalent roller. 


Because of symmetry of pressure profile about x = o, 
we have ^ = o at x = o. Using this condition in the 
integration of eqn. (3.28) we get# 

. -m, Vx )n Uj 2 n+l-q ( 3 , 30 ) 

dx In If ) h, h 

o i 

Integrating eqn* (3.30) with condition p = o at x = d# 
we obtain the expression for pressure. Denoting it by q 
we have 




A/ 


1 X 


(f )’^ h' 
o 


(3-31) 


The load component is given by 


«K,q = " 


(x) dx 


(3,32) 
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which on using eqn. ( 3 . 3l) simplifies to 

c 

/ 


h- 2m, v)n (1 ,q / - dx 

K, q 1 n h^ ^ (f 


(3.33) 


The time of squeezing from an initial film thickness 

q 

2hj^ to a subsequent film thickness 2 h 2 is obtained by 

substituting -V = ^ in eqn. (3,33) as 

In 

. _ 2n+l / ^"^l a/n ,1 ^q/n . ^xl^dh 

"K,q - - n- ( o ,,2n,l.q 


To study the effects of consistency ratio K and 
thermal factor q on loud as well as response time f 
we define the following quantities I 


"K 


dH, 




w,. 

^(a^) K, 

II 

i 

II 

0 ^ ^ 

Wq 

1 , q 

^(a^)l,< 

1 

1 

!^2 i 

t. 

i l/n 

1/q 

»2 ^ 


dH 


t = 

q 


b 

II 

^0 

^ (a ) K, q 
„ o_ . . 

^(a^)K,o 

K,0 

1 

0 

*K«q ^ 

/ 

""2 

l/n 

^(a2_)K,q 

^,0 

1 

_l/n 


/ 

^(aj^)K*o 


dH 


dH. 


(3,34) 


(3 . 35 ) 


(3.36) 


(3.37) 


(3-38) 


H. 
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where 


(a . )K, q 


1 

f 

o 


\n 


[x"+Vf %r !]« 

i 


2n-M— q 
^(a . ) 


i= 0, 1, 

(3.39) 


))= l-d-K”^/^) { ) )(2n+l)/nj (3.40) 

X = x/d, aj_ = a/h^, = d^/2Rh^ i = 1, 2, = h^/h^ 

1 H(a,) = % * <3.41) 

o 1 

From Eqns. (3.35) and (3.3 6) and tj^ are 
calculated to determine the effect of K . They are represented 
against K in Figs, 3.9 and 3.10 respectively- Due to ' 
consistency variation across the film thichness # for 
K <1 there is a decrease in the load ratio and response 
time for pseudoplastics relative to dilatants, and for 
K>1 these characteristics increase for pseudoplastics 
compared to dilatants. Hence ^we conclude that the change 
in consistency in the peripheral region alters the load 
and response time for pseudoplastics compared to dilatants. 
Similar behaviour is observed for q = 0 in the reference 

[is]. 

Figs. 3.11 and 3.12 depict the variation of load 
ratio and response time ratio with factors and 
for various values of n. For fixed value of h^ and h^/ 
and can be interpreted to represent the variation 
of peripheral layer thickness. Thus, we observe from 
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Fig. 3.11 that an increase in peripheral layer thickness 
increases the load ratio for K > l and decreases it for 
k <1» The effect of peripheral layer for K> 1 is to 
increase the load capacity more for pseudoplastics compared 
to dilatants. The trend is reversed for K< 1, Similar 
interpretations can be made for response time ratio tj^ 
vs. (Fig. 3, 12). 

Calculations of W and t from eqns. (3.37) and 

( 3 . 38 ) help to study effect of q on these characteristics, which 

are tablulated against q in Table 1. We observe that 

the decrease : in the load ratio and response time with 

increasing q is enhanced for pseudo plastics. 

Ta ble 1 


Variation of W and t with q for various values of n 
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3.5 JOURNAL BEARIKIGS 

The problem considered is the flow of a power law 
lubricant between two eccentric cylinders which approach 
each other with a relative squeeze velocity V, The velocity 
V is assumed to be constant both in magnitude and direction 
and symmetrical with respect to the boundaries of the 
bearing system. The clearance (difference in the radii) 
between the cylindrical surfaces is small compared to the 
radius of the inner cylinder. The configuration of the 
system is as shown in Fig, 3. 13* 

With the usual assumptions of lubrication theory, 
the governing eqns. of motion for a power law lubricant , 
in one dimensional form are given by 


dx ay 


ifo.au 

ay 3 y 


(3.42) 


+ .3 ^ = 0 

3 X ay 


(3.43) 


To determine velocity u from eqn.. (3.42) appropriate 

3 UL 

sign is to be attached to the velocity gradient • 

For this purpose let us assume that is zero at a film 

height y = h (x)^ then, in the region 0 ^ y / 

I'y^ 0 and in the region h < y < h, | ^ < 0. 


Considering the region 


^ < 0 , 
3 y - 


(3. 42) can be 



3 UnH , 

3y^ ^ 


written as 


(3.44) 
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Integrating egn. (3.4 4) -vd-th the boundary conditions 


= 0 at Y = h 
u = O at y = h 

the velocity profile is obtained as 


(3*45) 


u 


y 


dy h < y < h 
(3-46) 


3u 


Similarly^ considering the region ^0# we have 


u 




m 


dy O < y < h. 

(3-.47) 

The velocities are continuous ’ at y = h y hence from 
eqns . (3*46) and (3-47) we get 

■iW *** 

h y-h , , .. h -Y , , 

I (- )^-^ dy = / ( ) ' dy (3-48) 


m 


(3.4 9) 


h ■■■ o 

where the consistency m is given by (Fig, 3* 13) 

m = Km^(h/h^)^ 0 < y < a 

= m 2 (h/hj^)*^ a < y < h-a 

= Kmj_(h/h^)^ h-a < y < h. 

Using eqn, (3-49) in eqn. (3.48) we get h, = h/2 for a< h 
The volume flow flux of the lubricant is given by 


Q = / u dy 

o 


(3-50) 
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Using eqns. (3.46), (3.47) and (3.49) in eqn. (3.50) 

we get 


n - ■ » (il (2n+l)/n (2n+l-q)/n / \ / 1 dP'j 1/n 

° - 2n+l- ‘■1^ ■ \ h %) <- — gH) . 

( 3 . 51 ) 

where 

(g^) = l~(l-K"^/’^)[l-(l_2a/h) 

- { (1- ^)-(l- ^2n+l) /n (3,52) 

Integrating the eqn. of continuity (3.43) with the boundary 
conditions 

V = V at y = 0 

(3.53) 

= 0 at y = h 


we get 


•1^ = -V (3.54) 

9 X 

Using the expression for Q from eqn. (s."^!) in eqn. (3.54) 
we obtain 


d r 2n /l^ (2rH-l)/n , q,/n , (2n+l-q)/n/ ) r _ , 

dxl-2n+l’”4^ ^1 ^ dx -• 

(3*55) 

Substituting x = R© / where R is the radius of the inner 


cylinder, we have 
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i. (1^ t2n+l)/n ^q/n , {2n+l-q)/n/ w 1 dpv 1/nn 

R m 2n+V2^ ' h':^ h R ^ ] 

(3.56) 


= -V 


where 


h = c(l-ecos0) , h^= h(e = n) =c(l+e) 
V = -c || cos9 


(3-. 57) 


Integrating eqn. (3.56) and noting that ^ - 0 at 0 = O we 
have/ 


de 


SE r 22n+l(l jq(2n+l ^ ^ 1 )n siR“9 

1 hj^ 2n dt ^Qq) j^2n+l.-q 


. n. 


(3.5B) 


Since the pressure profile is SYmmetric about 0= O or 0= ^ / 
it is sufficient to consider the region 0 < 9 < ^ , where 
^ _< 0, to calculate the bearing characteristics. Integrating 
eqn. (3.58) with the full journal boundary condition viz. 
p(7i) = 0, the expression for pressure p can be obtained. 

Denoting p by p,, we have, 

r,.o2n+l/l sq/2n+l .sin^„9, , .p 

%,q = ^ ^ 2n 9^ (g^)n ^2ntl-q 

0 < 9 < Tt 

The load capacity Wj^ is given by 


(3.50) 


H 


= 2 


K.q o 




K/q 


COS0 R d0 


(3.60) 


which on using eqn. (3.59) beccmes 




TT „,.„n+l0 


sin 


K,q 


a 'dt' 


(g )" h^n+l-<3 


d9 


(3,61) 
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where 




(3.62) 


The Squeezing time ^ for the surfaces to approach from the 
initial concentric position e = 0 to a subsequent position. 


e = e^# say, is given by 


ae ]Vnae 

^K.rr ^ o m,2nJ.T>rr J 


O O (g )'"h 

o 


2n+l->q 


"k,, “1 

KJ 

(3. 63) 

The effects of consistency ratio K on load and response 
time can be studied by defining the quantities and 
tp, as follows I 

(3.64) 


W I 

W = K/g = 

K W, I, 

1/ q q 


^ = 


where 


^1 1/n 

t„ ^ I„ de 

K/ q_ = o q - 

t. 1/n , 

■ ■/q / d£ 

o 


(3 . 65 ) 


I = sin^"^,^ 9 . 

K,q o (G^)n(^.j^ 039 ) 2 n+l-q 


. de 


(3 #66) 


(G^) = l-(l-K“^/^)[l-(l-2a/H) 

- { (l-2a/H) -(l-2i/H) ^2n+l)/n j j 


(3.67) 


H = h/c = 1-ecose , a = a/c 


( 3 . 68 ) 
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Similarly/ the effect of thermal factor q can be studied 
Qn load and squeeze time by defining the quantities 


w 


q 


t = 


= (i_)q ^ 

^K/O «1 "K/O 






1 T-l/n 


f I 


jC{ _ 1 ^q/n o_ 


ds 


^,o 


'H 


^1 -1/n , 
o 


(3.69) 


(3-70) 


where is defined in eqn. (3.66) and H = h, /c = 1+e 

q ^ 1 1- 


The quantities ^ ' ^q ^q ‘^.efined in eqns. 
( 3 . 64 )/ ( 3 . 65 ), (3.69) and (3.70) are evaluated for 

various values of n and depicted in Figs. 3.14-3.20. The 
effect of K is to increase the load capacity for all n 
and this increase is enhanced for pseudoplastics for 
K > 1 (Fig, 3 '. 14 ) . There is negligible difference in load 
ratios for Newtonian and dilatant lubricants as K(< l) 
increases ,* these ratios are represented by a single curve. 
A trend similar to that of load ratio is observed in the 
case of response time ratio vs, consistency ratio (Pig, 3-15). 


From Fig. 3-16 we observe that for full journal 
bearings/ due to the effect of q,the load capacity decreases 
for all n and the decrease is enhanced for dilatants 
whereas the decrease in the response time is enhanced 
for pseudoplastics as evident frcsn Fig. 3 . 17, In the 
case of half journal the effect of q on and tg is 
represented in Table 2 , It is seen that the decrease in 
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Fig. 3.16 Effect of thermal factor q on load ratio Wq for 
journal bearing . 
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the load capacity and response time is enhanced for 
ps eado plastics . 
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Table 2 

Variation of and t^. with q for half journal bearing 


for 

various 

values of 

n 






W 

q 


1 

q 

' ' - - 

■ 

q 

n = 0,5' 

n = 1,0 

n = 1,5 

n = 0.5 

n = 1.0 

n = 1.1 

0 

1.0000 

1.0000 

1.0000 

1,0000 

1.0000 

1.0000 

.2 

.9169 

.9178 

.9184 

.9653 

,9837 

,9896 

.4 

.8411 

.8426 

■i.8437 

.9325 

.9679 

.9794 

•6 

.7718 

.773 9 

.7754 

.9013 

.9525 

.9694 





... 


» t. — - - 


Fig. 3. 18 represents the variation of ^ with 
for full journal bearing. Increase in implies a 
decrease in the film thickness'* The load capacity 
decreases for K< 1 and increases for K > 1 for all n, 
however, the difference in the ratio disappears for all 
n as the film thickness sufficiently beccmes smaller. For 
K < 1, the decrease in the load capacity is enhanced for 
pseudoplastic behavio'ur of the lubricant- Similar trend 
is observed in Fig. 3.19 for ■t^^ vs, except that the 
Increase or decrease of response time due to the decrease 
of film thickness is enhanced for pseudoplastics even for 
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small values of tliickness* As the film thickness decreases/ 
the effect of q enhances the decrease in the load capacity 
for dilatants and enhances the decrease in the response 
time for pseudoplastics CPig.3-20). 

3.6 SPHERICAL BEARINGS 

Consider squeezing of a power law lubricant between 
two eccentric spherical surfaces of radii r and R {r>R)/ 
approaching each other with a relative squeeze velocity 
V. The velocity is assumed to be constant in magnitude and 
in direction and symmetrically placed with respect to the 
boundaries of the syston (Fig, 3. 21). The film thickness 
is given by h = c(l-ecosO) where c = r-R and e = e/Cr-R), 
e being eccentricity ratio* Following a precedure similar 
to that adopted in the previous section/ one can obtain the 
velocity , u of the lubricant as 

u = (- ^ dy <0/ h/2 < y< h 

dx y m 9 y 

(3.71) 

= (_ ^)Vn ^ (Vl-j)!/!! dy > 0 / 0 < y <h/2 

dx o^ -y — — 

(3.72) 

For a sphere the amount of lubricant passing through a 

conical element of the surface is given by 
■ h 

Q = / 271 R sine u dy 

o 


(3-73) 
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1/n 


Expressing eqns, (3.71) and (3.72) in polar coordinates 
and substituting in eqn. (3.73) we get 

Q = 2TlR( -|P-^-){i)<2n+lVn h‘2n+l-qVn(g^,3i„9( 1 

( 3 . 74 ) 

where 

(g )=l-(l-K"^/^)[l-(l-% (2n+l)/n_ 2n+l£ (2n+l)/njn 

o ^ h n+1 h h J 


(3*75') 

hj|_ being the inlet film thickness measured at © = ”/2 
The displacement of lubricant at any point in the seat is 
given by [17] 

2 ^ 

Q = 7T VR sine (3,76) 

Prom eqns. (3.74 ) and (3.76) we get 


dp _ p /2n+l VRsn ~2n+l 

as - ' 2 ,; ■ 2 


sln^e 


fl nQ 


where V 



(3.77) 


Since the pressure profile is symmetric about ©= 0, we 
shall consider the region 0< 9 < ^/2 where ^ < O- 
Integrating eqn. (3.77) with boundary condition pi'^/2) = 0 
we obtain the pressure. Denoting it by p^,^ ^ we have 


( 0 ) 

%,q 


m^R <2 


Rc d^ 2n+ljn 


dt 2n 


22n+l ^1 


0 


71 /2 . n 

yq ' sin 


de 


0 < 0 < 71 /2 


(3.78) 
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The load capacity for honisphere is given by 

# M 




71 /2 / , 

K.q = 271 R ^ sinQ cose d9 

O 


which on using eqn* (3.78) yields 

K,q 2^h^^ ^ jnj^2n+l-q 


% 


(3.79) 


i 


(3.80) 


where 


2 1 2 2n 


(3 .81) 


The squeezing time for the surfaces to approach from the 
initial concentric position e = 0 to a subsequent position 
e = Ep, say , is obtained as 

71/2 

1 / T rT/« r / sin^"^^ 0 

^ = ( 2 )l/n(L.)Vn f [ •* / ‘ ,nj2nVl- 

^ ^ ^ o (g^) h^ 


sin— 9,_: 

q 


(3.82) 


As in the case of journal bearings^ we may define the 
following quantities to study the effects of K and q on 
the squeeze film lubrication of spherical bearings . 


w = 




^f/q - 

f^/q 

""l.q 

^l.q 

> 1/n 


^ Jt, 

de 

.P....K,q 


jl/n 

t A,q 

de 


(3.83) 


(3 »84) 
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and 


W. 




where 


J. 


w 



0 




^1 



/ 

o 

jl/n 

■^,o 

^1 



/ 

K,q 


o 


2 . n+2 

0 

2n+i-q 


_ . ^sxn 

o 

dO 


de 


(3.85") 


( 3 . 86 ) 


(3 - 99) 




( 3 - 100 ) 


i = a/c, H = h/c = l-£cos0 , 


3 .7 CONCLUSIONS 

In this Chapter, the characteristics of various 
squeeze film bearings with power law lubricant has been 
investigated taking into account of consistency variation. 

The effect of high consistency peripheral layer is to 
increase the load capacity and response time. The increase 
is more for pseudoplastics compared to dilutants. For 
the case of low consistency peripheral layer fluid film 
the achvex properties are reversed. As the periph 
layer thickness of high consistency increases the load 
and response time increase as has been shown in 
of roller bearing- 



The effect of the themal factor is to decrease 
the load capacity and response time. The decrease is 
more pronounced for pseudoplastics conpared to dilatants. 
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NOMENCLATURE 

peripheral layer thickness 

non-dimensional quantity corresponding to a 
radial clearance 

length of parallel plates/film stretch 
in roller bearings 
orh film thickness 

initial and subsequent film thicknesses 

minimum film thicknesses 

consistency ratio 

consistency index 

flow behaviour index 

hydrodynamic pressure 

thermal factor 

polar coordinate 

radius of the circular plate, sphere, jotirnal/ 

equivalent radius of the roller 

time 

response time 
response time ratios 
normal velocity 
load capacity 
load ratios 
eccentricity 

final eccentricity position 



CHAPTER 4 

EFFECTS OF CONSISTENCY VARIATICM FOR LIGHTLY 
LOADED ROLLERS IN COMBINED ROLLING 
AND NORMAL MOTION 

4.1 INTRODUCTION 

In Chapter 2 , we studied the lubrication of rollers 
under pure rolling. But in the actual bearing performance, 
usually there are transient and periodic forces or displacements 
imposed on the operating lubricant film and the rollers 
are subjected to both rolling/sliding and normal motion. 

Sasaki [ 1 ] was one of the earliest investigators to study 
the dynamic behaviour of roller bearings, Dowson et-al [2,3] 
made a theoretical study followed by an experimental verificat- 
ion of combined rolling /sliding and normal motion of rollers 
using a Newtonian fluid. In the case of non - Newtonian 
fluids, most papers which have appeared hitherto deal with 
either squeezing or rolling/sliding motion [4-9]. However, 
recently Sinha et.al [lO] theoretically analysed the lubrication 
of rollers in combined rolling and normal motion by character- 
izing the non— Newtonian behaviour of the operating lubricant 
by power law model » They did not give attention to variation 
of consistency in the lubricated region. It was pointed 
out by Qvale and Wiltshire [ 11 ] that lubricants would 
have a viscosity near the bearing surfaces differait from 
that of the bulk fluid owing to reaction of additives 
and surfactants with the bearing suf faces. As fluid 
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surface response to thermal effects is rather slow with 
time, thermal equilibrium, perhaps, be improbable with 
highly transient conditions in combined rolling and normal 
motion^ The principal objective of this Chapter is to 
study the qualitative behaviour of combined rolling and 
normal motion under the assumption of thermal equilibrium- 

Keeping this in view, in this Chapter, the effects 
of consistency variation due to aforesaid factors are 
considered across as well as along the film thickness 
in the case of lightly loaded rollers with a power law 
lubricant in combined rolling and normal motion considering 
cavitation- The consistency variation is accounted through 
the model suggested in Chapter 2 (eqn. {2-14) ) - 

^^2 BASIC EQUATIONS 

Consider the synmetrical flow of a power law lubricant 
between two identical rollers of radius r moving with a 
velocity U and a normal velocity V(Fig.4*-l)* The pressure 
and velocity profiles are depicted in Fig. 4-2. Takxng 
into account of the consistency variation of the lubricant 
film, the' governing equations to determine pressure in 
the fluid, region are given by eqns. (2-17) and (2.18). 




y 

A 



5^a .0 

dx 


Fig. 4.2 Lubrication of two rollers with power 
law lubricant. 
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dx*- 


-E . (h ^2n+l)/n 
2n+l 2^ 


^(2n+l-q)/n^^ -i 

o m, J 


U dh 
2 dx 


+ V 


-OO < X < -3^ (4. 1) 



R (1) (2n+l)/n J,q/n 
2n+l 2 


^(2n+l-q)/n 

o - 



“^*1 ^ i^c 

where P 2 _ and P 2 denote the pressure p in the regicais 

defined in eqns. (4»l) and (4.2) respectively, -x * and x^ 

are the points of maximum pressure and of cavitation, 

hj_ is the inlet film thickness at which the thermal factor q is 

■determined and the function (f^) is given by eqn. (2-16)! 

(f^) = 1-(1-K“^/^) { l-(l-2a/h) (4.3) 

where K is the consistency ratio and a is the peripheral 
layer thickness. 

The film thickness h is given by 

h = h^ + xV (2R) (4.4) 

where h^ is the minimum film thickness and R = r/2 is the 
radius of the equivalent cylinder. 


To determine pressure frcm eqns. 
we use the following boundary conditions 


(4-1) and (4.2) 

(aqns. (2.2l)— (2.24) 


Pi = 


0 at X = - oo 


(4.5) 
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dp, dp. 


dx dx 


Pi = P' 


= 0 at X = - 


X 


at X = -X 


dPc 


P2 = dx ■ = ° = 


X, 


{4»6) 

(4.7) 

(4.8) 


Integrating eqns- (4.1) and (4.2) with condition (4.6), 
we obtain the integrated form of the Reynolds equation as. 


[ ^ (h-l^ ) +V ( x+x*) 


dx 


(2nti)n 2^"+^ m, (^)‘J 


n 


I'h,"' jn j^^n+l-q 


-OO < X < - 3 ^ 


(4.9) 


♦ \nn 


, [ ^(h*-h) -V ( x+x* ) ] 

T2 = .(2ntl)n 22ntl a )q 

.X n 1 h. 


dx 


h (f 

o 


-x‘^< x < x^ 


where h - h at x = -x . 


(4.10) 


Integrating again eqns. (4.9) and ( 4 .I 0 ) using conditions 
( 4 . 5 ), ( 4 - 7 ) and ( 4 . 8 ), we obtain the expressions determing 
pressure t 


X rU 


Pt(x) = 22^+1 m^(|-)^ f ^ 

1 -00 j 


[^(h-.h *) +v(x+x *) ] 


n 


n , 2n+l-q 


• dx 


-00 < X < -X 


( 4 - 11 ) 
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P2 Cx) 


/2n+l\n 
' n' 



C j (h-h *) +V ( x+x ) 
(i ' h2“+i-‘* 


dx 


t I (h*-h)-v(xfx*) f , 

~ ' - - . - ax} < 

-X* (f^)^ j^2n+l-q 


X < X 

— c 


(4.12) 


Using condition (4.8) in eqn. (4*12) we obtain 
-X Xc 

f I3 (n,K,q,Q)dX = / (n^K, q, Q) dX (4*13) 

-X* 

where 

I3 (n,K, q, Q) = [ (H-H*)+2Q(x+X* ) (4^ I4 ) 

I^(n,K^q^Q) = [ (H*-H)-2Q (X+X*) (4*15) 

(Fq)= 1-(1-K"^/^) { l-(l«2a/H) (4*16) 

X = x/'^2RhQ., X*= x^Y2RhQ/ X^= x^/f2Rh^, H=h/h^= l+X^/ 

H* = =\/(2R/hQ>;^,i = a/h^ (4 ‘I?) 


A linear relation between X and X^ can be obtained frcxn 
eqns* (4.8), (4.IO) and (4.17) as 

X = X*-2Q (4. 18) 

c 

where Q is the dimensionless normal velocity. 


The (normal) load w(n,K,q,Q) is given by 

r^c ^ dp j 

W(n,K,q,Q) = / pdx=-J x ^ dx 

-.00 

which on using eqns. (4.9) and (4 .10) becomes 



-X 


W(n,K,q,Q) = 22n+l j. 


/ xi^(h-h )+V(x+x ) f* 

(f Jn h2n+l-.q 
o 


dx*- 


x{|(h*-h)-v(x+x^!'‘ 

,* (f 

-X O 


dx] 


(4.20) 


To study the effect of consistency variation across the film 

thickness on load, we define the load ratio as follows I 

K. 

w(n, 1, q, Q) 


-X 


X, 


/ XI 3 (n,K,q,Q)dX - / XI^ (n,K, q, Q)dX 

“•X 


V 

/“^KL XI 3 (n, l,q,Q)dX-/^cJaxi (n,l,q,Q)dX 


-X 


* 

K1 


(4.21) 


The case K = 1 yields the results corresponding to no consistency 
variation across the film thickness*- that is^ the results 
pertaining to consistency variation along the film thickness 
only for specified values of n, q/Q and for the same instantaneous 
minimum film thickness as that of the case with the specified 
consistency ratio. The points -Xj^ 2 _ snd correspond to' 

points of maximum pressure and cavitation respec-fcively, for the case 

K = 1. 

The effect of consistency variation along the film thickness 

on load can be studied by defining the quantities as, 

^ _ W(n^Ky^q,Q)^ 

q w(n,K,0, q) 


-X 


f XI_ (n,K, q, Q)dX.- /* XI^ (n,K, q, Q)dX 


^X - 


(4 .22) 


^1 /"^qo XIv(n,K,0,Q)dX -/cqo'xi^(n,K,0,Q)dX 


-X 


qo 
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^ ’th 

where ~^qo ^cqp points of maximum pressure and 

cavitation in the case of q = 0 with the specified values 
of n/K, Q and with the same instantaneous minimum film 
thickness and H^=h/h^» 

Similarly, the effect of squeeze velocity Q on load is 
studied by defining the quantities as 


W (n , K, g, Q) 
“ w(n,k/q#0) 


—X 

f . Xl3(n,K,q,Q)dX -/ XI^ Cn,K, q,Q)dX 

-OO _J* 


( 4 . 23 ) 


/ 


-X 


Qo 


XI. (n, K, q, 0 ) dX- / XI . (n, K, q, 0)dX 


-X 


Qo 


where ~X^ and X are points of maximum pressure and 
Qo CQO 

cavitation in the case of pure rolling (Q=0) with the same 
instantaneous minimum film, thickness. 


The frictional force f 2 ^n,K, q>Q) on the surface 
h/2 is given by 

I. 

fo(n,K,q,Q) =/ ° t; ^ ^ ^x 
-00 ; y=h/2 


(4-24) 


where 


5u ^n-l 


9y 


^ ^ i ay 

Using the expressions for velocity u in different regions 
obtained in egns, (2,5) and (2,6) we get. 


(4.25) 


-X . dp^ dp. 


= -[/ t 




2 dx 


-X 


(4.26) 
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substituting the expressions for and from eqns.U .9) ana(4 .10) 
In eqn. (4.26) we study the effects of consistency variation 
across and along the film thickness and normal velocity 
Q on frictional forces by defining the quantiUes 

2K' 2q 

and respectively as follows ! 


2K 


f 2 ^n/ q) 

■f^ri/l^qVQ) 

X^l3 (n,K,q,Q)(aX-j^° X^I. (n^K^ q^Q) dX 


-X 


X^I (n, l,q,Q)dX-/ (n,l,q,Q)dX 


(4.27) 


-X 


Ko 


2q 


f 2 (n, K, qt Q) 
■2 


■ 2 

f.n (n,K,0,Q) 


1 

H? -X* 


X^I^ (n^K, q, Q) dX - X^I^ Cn, K, q, Q) dX 


-X 


(4.28) 


1 J qo x^l^in,K,0,Q)dX-f '^‘*^X^I^(n,K,0,Q)dX 


-X 


f2(n,K,q, q) 
2Q ~ f^^n, K, q, 6J 
-X 


10 


/ X^I- (n,K^q/Q)dX -/ # X^I * (n,K, q, Q) dX 


(4.29) 


f X^I. (n,K,q,0)dX -/ X^I. (n,K, q,0)dX 


-X 


Qo 


where the expressions l3(n/K, q, Q) and I^(n/K, q, Q) are defined 
in eqns . (4.14) and (4«15) respectively. 

/ , 
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4.3 RESULTS AND DISCUSSIONS 

Eqn. ( 4 .I 3 ) is numeJricHlly evaluated to determine 
the location of cavitation point for various values of 
Q and n. Positive values of Q represent the normal 
separation of rollers and negative values denote the 
squeeze velocities. However, in this study, the squeeze 
velocity alone is considered as is the case in most 
normal bearing performances. 

Pig3‘4«3 depicts the effect of squeeze velocity on 
the location of cavitation point for various values of n. 
As the squeeze velocity increases, the cavituation point 
moves av/ay from the minimum film thickness. The rate 
of shifting of cavitation point with respect to squeeze 
velocity is almost the same for all n except that for 
pseudo pi as tics it lies further away from the converging 
gap than that for dilatants. 

In Fig. 4.4-4.11 are studied the variation of load 
ratio and frictional drag ratio parameters for different 
n. As in the case of pure rolling, we observe that for 
specified squeeze velocity, an increase in a results in 
an increase in the load capacity for the case of high 
consistency peripheral layer fluid film (Fig. 4.4). 

A similar interpretation can be given to frictional drag 

ratio parameter ^2a 


X 
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In Fig. 4.6 is studied the variation of W,. with K 
for various values of n. Again we note that the effective 
increase in the high ccwsistency peripheral layer is to 
increase the load capacity for all n. The increase due 
to the effect of this peripheral layer is to give more 
load, for dilatants , however, in the case of low consistency 
peripheral layer more load is pronounced for pseudo pi as tics. 
Similar behaviour is observed for frictional drag ratio 
parameter (Fig, 4.7). 

In Fig, 4.8 is plotted the load ratio against squeeze 
velocity Q for different values of n. The effect of the squeeze 
velocity is to increase the pressure generation in the 
lubrication regime* Due to this effect the load capacity 
increases for all n. Also, it is reconfirmed that in the 
case of combined rolling and squeezing the load is more 
compared to pure rolling'. The increase in the load capacity 
due to squeeze velocity is more pronounced for large values 
of Squeeze velocity. The increase in the load capacity 
is more for dilatants. A similar interpretations can be 
given to frictional ratio parameter ^2Q Q(Fig.4.9). 

As in the case of pure rolling, we observe that the 
effect of the thermal factor q in combined rolling and 
squeezing is to reduce the load capacity and frictional 
forces (Figs. 4 . 10 and 4.11) . The reducticai is more pronounced 
for dilatants compared to pseudoplastic behaviour of the 


lubricant. 



4.4 CONCLUSIONS 

The effect of the squeeze velocity is to increase 
the load capacity and frictional drag for all values of 
the flow behaviour index n. This is because an increase 
in squeeze velocity means more pressure generation 
in the lubricated region. The increase is more for dilatants 
than that for pseudoplastics, as in the case of pure rolling. 
As the squeeze velocity increases, the cavitation point 
shifts away fron the minimum film thickness. 
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Fig. 4.9 Frictional drag ratio parameter F 2 Q vs. 

squeeze velocity for various values of 
flow behaviour index n. 





10 Load ratio parameter Wq vs. thermal factor 
q for various values of flow behaviour index n 
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nomenclature 


a peripheral layer thickness 

a non-dimensional quantity corresponding to a 

±2 frictional force at y = h/2 


^2K'^2q'^2Q Frictional force ratio parameters 


h 


h 


o 


h. 


H-.,, 


K 

m,mj^ 

n 

P 

q 

Q 

r 

R 

U/ V 

U 

V 


film thickness 
minimum film thickness 
inlet film thickness 

non-dimensional quantities corresponding to 

h and hj^ respectively 

consistency ratio 

consistency indices 

flow behaviour index 

pressure 

thermal factor 

non-dimensional quantity corresponding to V 

radius of the roller 

radius of the equivalent roller 

velocities along the coordinate axes 

rolling velocity 

normal velocity 

non-dimensional load ratio paraneters 



CHAPTER 5 


LUBR^CMvf^ VARIATION OF BOWER LAW 

LUBRICAivTTS FOR riqid ROLLERS WITH ROUGH SURFACES 


5,1 INTRODUCTION 

In earlier chapters while studying the effects of 
consistency variation of power law lubricants on bearing 
characteristics, an implicit assumption was that the 
bearing surfaces were perfectly smooth- But, it has been 
recognized that owing to machining limitations this 
assumption is rather xinrealistic. The topography of bearing 
surfaces is generally rough. It consists of a net work 
of randomly distributed micro-metallic structures forming 
^ hills and valleys' running on the surfaces. The study 
of these structures and their participation mechanism 
in the lubrication process is vital in determining the 
bearing characteristics, particularly when the dimensions 
of these structures are comparable with the thickness of 
the operating lubricant film. 

Surface roughness in lubrication problans is studied 
by following deterministic and stochastic approaches. 

Burton [ 1 ] and Wildmann [ 2 ] studied the effects of surface- 
roughness through deterministic approach by postulating 
a sine (cosine) wave or a series of sine (cosine) waves 
for the film thickness. Recently, a new deterministic theory 
has been proposed by Shukla [ 3 ] applicable when the mean 
height of the roughness structure is of the same 
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order as that of the minimum nominal (compliant) film 
thickness^ This theory for hydrodynamic or mixed lubrication 
proposes a division of lubricant flow region into roughness 
interaction zone along the rough surface and purely 
hydrodynamic zone along the smooth surface. In the stochastic 
approach, the film thickness is assumed to be a random 
quantity varying with tha roughness pattern. Christensen 
and Tender [ 4 ] presented a stochastic theory for hydrodynamic 
lubrication of rough surfaces with roughness lay running 
along and transverse to the direction of motion, by making 
two heuristic assumptions about the pressure gradient and 
fluid flow rate. It was concluded that the type of roughness 
along the direction of motion results in a slight decrease 
in load carrying capacity and an increase in frictional force 
thereby causing a significant increase in coefficient of 
friction. The effect of the roughness structure transverse 
to the direction of motion is, however, to improve the bearing 
characteristics. Christensen [ 5 ] extended the stochastic 
theory to mixed lubrication regime. Raj and Sinha [ 6 ] canployed 
Christensen' s model to analyse the effect of roughness in 
short journal bearin<^ assuming that the standard deviation 
of the distribution describing asperity heights and minimum 
nominal film thickness were of the same order. Rhow and 
Elrod [ 7 ^ considered the spatial characteristics of roughness 
with respect to amplitude and wave length claiming some advantage 
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over stochastic analysis of Christensen and Tender 'in both rigor 
and generality though not in brevity' . An average Reynolds equat- 
ion for three dimensional isotropic and nonisotropic roughness 
structures has been derived by Patir and Cheng [ 8 ] by defining 
pressure flow factors obtained through flow simulation. Recently, 
Tender [ 9 ] developed a roughness theory based on numerical 
evaluetion of average flows and shear stresses. 

The study of surface roughness together with the rheological 
considerations of the lubricant has been undertaken by some research 
workers. For e.g.,, Shukla and Kijmar [ lo ] and Kumar and 
Sachidcinanda 11 ] considered effects of viscosity variation and 
surface roughness. Thermal effects in the lubrication of rough 
surfaces were studied by Dyson [12] and Fowles [13]. 

With the development of theoretical models to account for 
surface roughness in lubrication problsn, a good deal of 
experimental work on surface roughness has been reported’. 

Tsukada and Sasajima [ 14 ] developed a measuring systen to 
investigate the three dimensional characteristics of 
asperities.. An interf erometic study of rough surfaces under 
ehd lubrication was conducted by Jackson and Cameron [ 15 ]. 

An extensive review of experimental report of various techniques 
for the study of separation of surfaces, the real area of contact. 
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the number of contact spots, the spatial distribution of 
the contact spots, the distribution of their sizes and the 
relation of all these to roughness and to the normal load 
was presented by Woo and Thomas [ 16 ]. 

Despite the development of various models to include 
the effects of surface roughness in lubrication problans / 
a critical analysis reveals that a most practical and general 
model is yet to emerge [ 17- ]. In reviewing the Christensen's 
model /Sun [ 18 ] questioned the validity of the heuristic 
assumptions made about pressure gradient and flow rate and 
pointed out the absence of roughness correlation terms in 
this model. Chen and Sun [ 17 ] ranarked that Elrod's model 
[ 19 ] suffered from the restriction that the roughness spacing was 
small. The three dimensional roughness theory of Patir 
and Cheng [ 8 ] has the potentiality to describe the arbitrary 
roughness arrangement , However, the flow factors determined 
in this model are subject to alterations by the choice of 
numerical representation of contact, shape, size and boundary 
conditions of the model bearing, as has been danonstrated by 
Teale and Lebeck [ 20 ]. Despite the complexity or stortcomings 
of various models for surface roughness, the study of the 
rheological properties of the lubricant on rough bearings was 
attempted. Shukla and Kumar £ 21 ] considered viscosity 
variation and surface roughness in the analysis of slider 
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bearing lubrication using Newtonian fluid. Prakash et.al. [22] 
studied, the roughness effects in slider bearing using ( 3 . 
micropolar fluid and pointed out that the effects like 
enhancement of viscosity in thin film zones were rheological 
effects and might not be attributed to the surface roughness. 

In this Chapter, we study the effects of surface 
roughness and consistency variation for rough rollers using 
a power law lubricant. Such a variation in consistency may 
be warranted by the reaction of additives and surfactants 
[ 23/24 ] .C onsistency variation is accounted through the 
model suggested in Chapter 2 /while surface roughness is 
considered through Christensen's stochastic theory £ 25 ]• 

5 \ )T0CHASTIC F0RI4o OF REYNOLDS EQUATION 

The problem considered is that of the lubrication of 
two identical rollers each of radius r, lubricated 
with a power lav; fluid(Fig.5,l) . U is the rolling velocity 
and V/2 is the normal velocity of the rollers. The rollers 
are separated by a film thickness h which is considered 
to be a stochastic variable depending on the roughness structure 
of the bearing surface • FurtheU/we assume the conditions for 
applicability of Reynolds equation for treating rough surfaces [17 

i. e./ 

(1) the roughness spacing is large catipared with the 

film thickness/ and i 





il 




A schematic diagram of a rou^ 
roller bearing. 
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( 2 ) the roughness height is small compared with 
the film thickness so that no recirculation 
region forms in the film. 

As pointed out earlier, the bearing surfaces are 
generally rough. To take this into account, the film 
thickness is to be viewed in the light of surface 
asperities in the liibricated region. Stochastically, 
the film thickness is expressed as [ 25" ] 

h = h^(x, 2 , t) + hg(x,z, ( 5 . 1 ) 

where h^ denotes the nominal part of the film geanetry 
and h is the roughness part measured fran the ncmiinal 

o 

level. ^ appearing in h is a random variable describing 
the roughness arrangement. Assigning a particular value 
to t means selecting a particular arrangannent from the 
supply of possible roughness arranganents. For a given 
value of the surface roughness component h^ is a 
deterministic function of the space variables. We 
consider for our study one dimensional roughness pattern 
formed in the lubricating machinery, particularly after 
they have received some run-in and wear [ 6 , 24]* Keeping 
note of this, the film thickness in the case of rollers .can 
be written as [ 8 ] 


h = h^ + 6 l + 62 


(5.2) 
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where 63 follow Gaussian distribution function mth 

mean zero and ■ combined standard deviation 0 . 


For many engineering surfaces the use of Gaussian 
function to describe the roughness profile at the vertical 


scale is valid at least upto three standard deviations [[ 26 j . 
For mathematical simplicity but with a fair degree of accuracy^ 
a polynomial function which closely approximates Gaussian 
function [ 27 ] is employed in the und'ergoing analysis T 


f (h ) = ((35/32) /b"^ )(b^-h 2)3 ^ < 


b 


= O 


(5.3) 


otherwise 


where the relationship between the roughness amplitude b 
measured from nominal level and stardard deviation a is 
given by b = + 30. 

In the following section/ we shall derive stochastic 
form of Reynolds equation considering cavitation, applicable 
to rough rolling surf aces, using Christensoa^ s approach [ 24 ]• 
For this purpose, the expected values of various quantities 
may be defined as follows * 

e( ) = /°° ( ) f(hg)dhg {5-.4) 

00 

where f(h 3 ) is the probability density function of the 

s‘ 


random variable h 
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If the function is assumed to be symmetric, th^ 

^ 2r+l 

f =0 r = 0,1,2,.. (5.5) 

This assumption is valid since the mean of h over the 

s 

bearing surface is zero and can be taken as the definition 
of nominal film thickness [ 24 ]. In particular, for the 
ergodic (stationary) stochastic process, we have 


E (h) = 

With the help of eqns. (2.10) and (2.11) arxi Fig. 5.1 
the basic equations determining pressure are written as 


where 


^ {(f) (?^i) Vn , ^_V U ^ -co< X < -X* 

dx dx ^ 2 2 dx ^ ^ 

(5j6 ) 

^ 5 (f) (- Vn , .(.V U ai) ^ 

dx ^ dx ^ 2 ^ 2 dx ^ ^ 

(5-.7,) 

h/2 T , 

' - (5^.8 ) 


(f ) = / y(y/m) ' ^ dy 


o 


and -X* and x^ are the points of maximum pressure and 
cavitation respectively ^p^ and P 2 denote the pressure p in 
the regions defined in eqns. (5£) and (5.7) where pressure 
gradients are positive and negative respectively. 
Considering consistency variation, the consistency m is 
given by 

m = mj(h/hj)*^ 0 < y < i]j/2 - ^-p 
= Kmi(h/h^)'^ A/2-^p < y 1 


(5. 9) 


I 
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where is the apparent peripheral layer measured from 
the nominal level* h^ is inlet film thickness where q is measured. 
Substituting the expression for m from eqn. (5.9) in eqn. (5.8) 
we obtain 

(f) = (n/( 2 n+l)) ( 1 / 2 ) ( i /m 3 _) ^ (f^) ( 5 -. 10 ) 

where (f^.) = { (l-K”^/’^) (h^ j^( 2 n+l)/n 

(5.11) 

It is clear that for the smooth case/(f) in eqn. ( 5 . 10 )is 
reduced to that obtained in eqn. (2,15), 

(i) ONE DIMENSIONAL LONGITUDINAL' ROUGHNESS 

In this type of roughness, the lay of roughness structure 
consists of long narrow ridges and valleys running parallel 
to the rolling direction [ 24 ] (Fig. 5,2 (a) ) . Thus, the 
roughness components 62 _ and 62 independent of x, U and 
t (time) [ 28 ] . We note that the pressure gradiait ^ 
is a variable with zero or negligible variance and we have 
[ 25 ], [ 29 ] 

E(+ g)Vn 4 a|0>,)Vn (5.12) 

Now, taking expectation on both sides of eqns. (5'.6 and 
( 5 . 7 ) ,WG obtain the stochastic foiro of Reynolds equation 
applicable to rollers with longitudinal surface roughness . 



/ V 



' 1 


X 



iM-»-J 4A fc.^ ■„ , K. .■I^itvkwrwi. 


m4JL^. z 

7/7777 n7777777777f 


Fig. 5.2 (a) Longitudinal roughness 



Fig. 5.2 (b) Transverse roughness. 
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_£_r_2n _ ^ 

ax^ 2n+l 


(i, (2n+l)/„ ^ «"‘P£,l/n 

m, dx 


3 ^ 

= U - 


3h 


n 


3 x" ■'■ at 


«.<X) < ^ ^yr' 


(5.13.) 


d [ -_2n_ /l\(2n+l) /n , Q/n .n. i 

diC ^ - ax 


dE(p„) , , 

^ -flv- > ] 


3 h a h 


(5.14) 


3h 
-ja. 


where The nominal film thickness is given by 


^ 1=^0 + ^ /^ 2 R ) 


(5-.15) 


where is the nominal minimum film thickness and R is 
the radius of the equivalent roller. 


In this Chapter, we study the effect of roughness 
on rigid rolling surfaces only. For this purpose, we have 
to solve the following equations to obtain average pressure 
E(p) I 


^„2n_ (1) (2n+l)/n ^q/xi g/.. ^®^^l^,^l/n-| 

d;^2n+i ^1 dx ^ J 


U 


dh 

. n 

dx 


- oo < X < -X 


d r 2n ,1^ (2n+l)/n ,q/n L> ‘^^^^^tl/n-i ^ 


(5.16) 

dh„ 


-X*< X< x^ 


dx 
(5.17) 


We use the following boundary conditions 
E(p^) =0 at X = 


.. CO 


(5.18) 
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dE(p ) 

dE(p2) 


"dF'' ■ ^ 

'’(ix' "■ = 0 at x = -X* 

(5.19) 

E(p^) = 

at X = -X* 

(5.20 ) 


dE(p2) 


E(p2) = 

” dx"' = 0 at X = x^ 

(5.21) 


Integrating eqns. (5,1.6) and (5.17) with condition (5.19) 
we get 


dE(p^) 

~~dl' ■ 


(.2n+l)n 22n+l „ ^(tjq 

■in 1 


(h -h 
n n 


Ce([^))} 


n 


-•“ < X < -X* ( 5 . 22 ) 


dE(pT) o _l.t V, o j-T 1 (h -h )^ 

ta 2n 1 j lAfj.));” 

-x^ < X < x_ ( 5 •23) 

— — c 

* * 
where h„ = h„ at x = -x . 
n n 

Using eqns. (5.15) and (5'.2i) in eqn, (5'.23) we obtain 
* . 

Integrating again eqns. (5.22) and (5.23), with conditions 
(S.rs) , ( 5 - 20 ) and (5.21) we obtain the Average pressure 


given by 

E(pi) 


(?u+l)n 22n+l yH 
2n J- 


”1 


X 

/ 

—00 


C E^f 

r 


dx 


_ 00 < X < —3^ 


(5'.2f) 
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E(p,) = 




{E((f^))P 


# 

/II —II V 

.(.n n) 


h. -h _ ^ n 


/ ■ dx 1, 

_^* {E((f^))}^ 


dx. 


(5.25) 


Using condition (5.21) in egn. (5.25), we get the eguation 
to determine cavitation point ; 

-X 

/ ^ = / 1l2 (5.26) 

- CC 


where 


Ij^^(n,K,q,B) = (H^^# )’^/ { e((p^))}' 

I^2^n,K,q,B) = (H^-H^)""/ { EC(P^))}' 

35 B 2 ? ^ 

eCCF ))=- „-/ (B -Hs^)^(P^)dHs 
^ 3 2b"_b ^ 


(5.27) 

(5.28) 

(5.29) 


(F^) = [ (l-K”^'^^) (H^-2ap) ^2n+l)/n 

^.K~Vn(H +H (5-.30) 

K: = x/'^(2“'Rh^), x'4*/ fWi^h H = H = h/h^, = h._/h, 

9 * # -jf 2 — 

"n = V^o = 1+"^ ' t = 5> /*’o = 1+=^ ' 


S' o 


B = b/h 


o 


p o' 
(S'iSl) 


The load component Wj^^(n,K,q,b) per unit length 
in the x direction for the case of longitudinal roughn^s 
is given by 


-x* 2 '^E(p^) 

WL^(n,K,q,b) = / x 


X* 2 dE(p2) 

a, ax+ , X 4, dx 

—X / c o o ^ 
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which on using eqns. (5.22) and (5.23) becones 


W. 


bx 


9^.1 ^ 9 TT^ ■, ^ X (h -h 

( v >.'1 _ r^n+l\n ^2n+l U ,1 Nq? . . n. xi . , 

(n,K, q,b) = 2 m, ^ (tT") l ’ r> 

2n 1 R h^ "-coCECf^))}'' 


-2 


x*xh#-h)" 

: lu -a-., ax ] 


^ { ECf ))} 

o r 


(5.33) 


We obtain the load ccmponent ■'‘by (n,K,q^b) along 
Y' diroction is given by 

by 


-X* x{h -h'*^)^' 

2n+l,n ,2n+l „ n (l_jq n 


WT.„(n,K,q,b) = ^ ‘“1" 'h 


(5.34) 

To study the effects of longitudinal roughness and 

consistency variation on load, we calculate the load ratios 

Wr and defined as follows \ 

LX 


w, ^ (n,K, q, b) = 

JUX 


-X 


\x(n' 


Ww (n,K,q) 

Sx 




/ (n,K,q,3)dX - 2 / X * 5 / S) dX 

bl r> 


-1 


* 

_x* 4 

/ x2l„,(n,K,q)aX-2 S X^Igj (n,K,q) aX 

.-.oo ^ 


(5.35) 


dx 



(n,K, q,B)dX 


/ ^L1 



/ Xlgj^ (n,K, q) dX 


— 00 


(5.36) 


where Wg^(n,K,q) and Wgy(n,K,q) are given by (eqns. (2.36) 
and (2.38)) 


H 


S^(n,K,q) = (.^il-)^ 22^+1 


(h -if 


dx 


-2 


X 3 * 


dx] 


(f )" 


( 5 . 37 ) 


and 


Wg^(n,K,q) 


2n 1 h. 




-xf x(h -h*)^ 
jTS n. 

-CO (f^) h^ ^ 


(5.38) 


where -y^ and Xg are points of maximum pressure and 
cavitation, respectively, for smooth case and Igj^(n,K, q) and 
Ig^(n,K, q) are given by 


3 ^ = (H„-H*)7l(Fo)'' 

(5 .39) 

52 “ <4-Hn)"/ f teo)” ' 

(5-40) 




(P^)being given by eqn. (2.33). The cavitation point Xg in the 
smooth case can be obtained from eiyi- (2-3l)» 



(ii) ONE DIMEl-ISI'JNAL TRANSVERSE ROUGHNESS 

In this case, the roughness lay consists of long 
narrow ridges and valleys transverse to the direction 
of motion | 25 ] (Fig. 5 .2Xb)).This type of roughness 
pattern may be generated in bearings subject to oscillations 
perpendicular to the principal motion [ 6 ]. The film 
thickness in the case of transverse roughness is given 
by [ 27 ] 


h = b +6 t (x-Ut) +69 (x-Ut) • 
n 1 I 

a 3 ^1 

Noting that S^Cx-Ut) = -U 

3 X 3 X 3x 3x 


i 


1,2, 


( 5 . 41 ) 


(5-42) 

(S.43) 


egns . (5-6 ) and (5.7" ) become 


d r 2n 

"dF 2n+l 


fl. (2n+l)/n , q/n 
^2^ ' ^1 


(f, 




'"’at' 


.4 r„ 2 a_ (i) (2n+l)/n (_ U - = 

dx L 2n+l 2 ■ 1 r m^^ cix 

( 5 . 45 ) 


-«> < X < -X* (H.44) 

dp^ 

1 


3 t 


-X < X < x^ 


We denote the square bracketed term in eqn. (5-.44-) by M 
in accordance with reference |28j« 

Th 

„ . (1) (2"+l)/n hf" (f J (i- 

^ - 2n+l 2 ' 1 ^ ^1 
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On rearranging, we have 




(- 2n 

2n+l^ 


/I I 2n+l 




(5.47) 


In this case/M is a stochastic variable with zero or 
negligible variance. Hence/ we have [29] 

E(M+Uh^)'' = {E(M)+Uh^f = (M+Uh^)"^ (5.48) 

Taking expectation for terms in eqn. (5>44'') using eqns. (5-4 6) 
and (5.48) we get 


2n ,|,(2n+l)/n ,a/n {^j;. 1 ^,-1/n (1. 
zn+i. z J- (f dx -i 


ahn ah 

” ^ "Sx*' **" 9 t -“’<x< -X* (5.49) 

Similarly taking expectation/ eqn. (^.45) beccmes 


d r 2n ,'In (2n+l) /n 
dxLzn+l 2^ 




1 ''n 


1 


8h ah^ ^ . / ..X 

=-(u -fr + “3^^ ^ 1 ^c 

The average P'-eynolds equation for the case of transverse 
roughness under pure rolling conc3itions becOTies 


^1^ (2n+l)/n ( h { E ( 

dx (■ 2n+l 2 1 (f )'^ 

dh ^ 

_ u - -CO < X < -X 

dx — — 



l/n, 

dx"'^ 


(5.51) 
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” ‘35 <x:„ >.’^. 52 ) 

**■ •“ c* 

Following the same procedure as adopted in the case of 
longitudinal roughness, the expressions for H(pj_) and E(p2) 
can be obtained as 

e(pO - A(h_-it5“E(-l ^.)ax 


1 - oo n n {f 


-00^ X £ -x' 

-X 


(5.53) 


E 


(pj) = -(2|U)n U--* ia,[ , (Vl^F 

1 .. 


-00 


J* e(-^- j^)dx] ^ (5.54) 


.-i^- 




The point of cavitation X* is obtained fron the ecyi. 
•ifr 


-X 


X 


J . I^^(n,K,q,B)dX = 


-X 


(5.55) 


where Tiji2^(n/K, q, B) — (Hj^— I^) 


(5.56) 


1,^2 

The load components W™^(n,K,q,b) and W^(n,K,q,b) are 
given by 
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W^^(n,K, q,b) 


/2n-l-l'jn ~2n+l u’^/1 “^'^9 ;<■ t 

“'Tn-^ 2 ”1 / x2(h„.^)='- E(-l 

"1 -» ^ ^ (fp" 


* 


)dx 


X 


~2 / x^(h -h; E(—- dx 1 

o n n if^)^ ^ 


(5.58) 


and 


(n,K, q,b) = 


(2^,n 22n-.l 


X* * 


h. 


/ x(h^~h^)“ E(^-i--)dx 

(f 

~ CO 3;;; 


( 5 . 59 ) 

As in the caSa of longitudinal roughness, we define the 
load ratios and to study the effects of consistency 

variation and roughness for transverse case.' 


Wyy_ (n,K, q, B) 


(n, K, q, b) -Wg^(n. K, q) 
Wg^(n,K, q) 


—X 

/ X^I^^(n,K,q,B)dX-2 / X^I^2 

-00 o 


—X X 

J. ®X^l3^(n, K,q)dX -2 X^Ig^ (n,K,q)dX 
-00 ° , 


1 


(5. 60) 


Wyy = 


Wy^(n,K,q,b)-W^(n,K, q) 


Wgy(n,K,q) 


-X 


-!«■ 


j. XIy^(n,K,q,B)dX 
-CO 


_Xc 

f XIg^(n,K,q)dX 

-00 


--1 


(5. 61) 
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In the absence of consistency variation (i.e. K = l 
and q = O) the load ratio quantities for longitudinal 
and transverse roughness reduce to those in the 
reference [30] . 

5.3 RESULTS MTD DISCUSSION 

'.'■■Jhili.; dividing the operating lubricant film into central 

layer and peripheral layer with different consistencies, we 

defined in eqn. (5-9) apparent peripheral layer a . This 

layer is to be distinguished from the actual peripheral layer 

* a' which is the sum of the apparent peripheral layer and 

the roughness part of the film thic]cness from the mean level 

on each surface- Precisely t'a' is stochastic and 'a ' is 

P 

deterministic quantities. It is obvious that it is the 
actual peripheral layer which takes part in the lubrication 
process. 

In the calculation procedure, the av.erage pressure is 

obtained by using a symmetric function which approximates 

Gaussian function. This necessitates the roughness amplitude 

b/2 at each rolling surface to be less than or equal to 

(i-e.- a, > b/2). This enables the ^hills and valleys' to be 

P 

completely encompassed by the apparent peripheral layer. Since 
the rougl'iness heights are measured frexn the mean level, a^ may 
be measured from the mean level. In the theoretical analysis, to 
study the effect of consistency variation across the film 
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thickness, can be regarded as a paraneter. 

To study the effect of longitudinal roughness and consistency 
variation, the cavitation points X for the longitudinal case 
and for the smooth case are determined fror. eqns. (5-26) 
and (2 # 3*1 ) vand the load ratios and evaluated 

using cqns. (5.35) and (5.36) respectively. Similarly :X^ 
for transverse case can be determined from eqn- (5,55 ) and load 
ratios and can be evaluated using eqns (s,60) and 

(5-6l) respectively. Table 1 gives the picture of the effect 
of roughness parameter B on the load ratios and for 
various values of n for the cases K> 1 and K < 1, Table 2 
gives the picture of load ratios different 

values of B. It is seen from the Tables land 2 that the roughness 


effect is to decrease the load capacity in the longitudinal 
case and increase it in the transverse case for all n for any 
fixed a and K. The results obtained are in good qualitative 
agreement with those of Sinha et.al. [ 3® ] . In the case of 
longitudinal rouglafiosS/ for K < Ij, there is slight increase 
in load capacity as compared to the ca.3e of K > 1 for all n. 

It is evident that along the direction of motion the hydrodynamic 
action of the lubricant contributes to the load carrying 
capacity . In the longitudinal case, both for K > 1 and K< 1, 
the peripheral layer is attached to the bearing surface and the 

lubrication process depends mainly on the con sis tcaacy of the 

central layer. If the consistency in this layer is higher than 
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that of the peripheral layer there is an increase in the 
load capacity. Hence K < 1 brings in an increase in the 
load capacity. 

In the case of transverse roughness^ for K > 1, the 
high consistency fluid entra.pped between the asperities in domin- 
ant providing resistance to motion, which in turn 
increases the load capacity. This dcminance is supplefnented 
by the predominrant squeezing motion of the fluid in the 
transverse case and by the extensive length of the bearing 
in the axial direction. Hence we have the increase in 
load capacity for K >, 1 for transverse roughness. 

Considering no consistency variation across the film 
thictoiess for the smooth case, evaluation of the load ratios 
enable as to study the combined effect of peripheral consistency 
variations and surface roughness which is depicted in Fig, 5.3- 
It is observed that as K increases, the load capacity increases 
for all n. The trend is similar to that of the smooth case 
discussed in Chapter 2. However, for seme value of K lying bet- 
ween O.A and 0.7, this trend is violated i.e. there is slight 
decrease in the load ratio for n = 0.5 compared to n = 1.0 
for K < 1 . Thus we see that there exists a critical value of K 
(say , I<^) v<rhich may depend on the flow behaviour index n as 
well os the roughness parameter B below which the load ratio 
trend is different from that of the smooth case. The effect 
of the rheological anomalies could be complex for these values 
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^ pattern for flow bahaviour 

index can be predicted. This is elaborated in Table 3 
v;here the load ratio and arc tabulated against 

different n for fixed K(=0..4) . 

VJe observe an increase in the load ratio n decreases 

fran 0.3 to 0.7. It pronably reaches a maximum at some value 
of n lying between 0,7 anc5 1.1 ard then starts decreasing. 

This deviation in the load ratios can be attributed to the 
complex nature of the rheological anomalies in the presence 
of asperities. Similar behaviour is observed in the transverse 
case as well. 


Table 3 


The load ratios 


and W, 


flow behaviour index n 


n 

W 


.3 

.743 

.744 

.7 

.840 

-.840 

1.1 

.827 

.828 

1.5 

.798 

■i799 

1-9 

.770 

,771 



K=0.4 a^=0.3 B=0.2 q=0.5 

P 


In the case of transverse roughness, similar interpretations 
(as advanced for Pig.5»3) can be made for the load ratios 
and 
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In Fig.. 5.5 is depicted the effect of K on VL„ and W, ^ 

IjX IjY 

for various values of q. For K < l as q increases, the load 

ratios increases for all n and for K > l load capacity decreases 

for increasing q low values of q effect the load ratios W,. 

Lx 

and significantly for any K. Similar interpretations 

can be given for the case of transverse roughness (Fig. 5, 6 ) » 

5.4 CONCLUSION 

With the consideration of consistency variation across 
the film thickness^ the effect of the longitudinal roughness 
is to re-uce the load capacity for all n and transverse roughness 
is to increase it. Because of the interaction of the peripheral 
layer with the surface asperities, there exists a critical 
value of the consistency ratio above which the pattern of the 
flow behaviour index is similar to that of the consistency 
variation in the smooth case-. When the consistency is higher in the 
central layer, there is greater lodkl capacity for longitudinal 
roughness^ However when the consistency is higher in the peripheral 
region in the transverse case, there is greater protection 
against seizure and load capacity » With the given consistcHicy 
variation parameter across the film, the effect of the thermal 
factor is to increase the load ratio for a decrease of 
consistency in the peripheral region and decrease the load ratio 
for increasing consistency in the peripheral region for 


all values of n* 




Effect of consistency ratio K on load 
ratios Wlx various values 

of flow behaviour index n 





Fig 5-5 Effect of consistency ratio K on load 
ratios Wlx , Wly various values of 
thermal factor q 




K 


Fig 5-6 Effect of consisjency ratio K on load 
ratios Wrxand Wty for various values 

of thermal factor q 
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m/mj 


n 

P/ P1/P2 
q 

r 

R 

t 

U 


u, V 


■^Lx'% 

^Ty 

Wsx^"''^Sy 

\-X'^LY'^^TX'^''^TY 


apparent peripheral layo: 
roughness amplitude 
expectaticn of ( ) 

film thickness 

nominal part of film thickness 

roughness part of film thickness 

consistency ratio 

consistency indices 

flow behaviour index 

hydrodynamic pressures 

thermal factor 

radius of the roller 

radius of the equivalent roller 

time 

rolling velocity 

velocities along the coordinate axes 
loads in the longitudinal roughness 
loads in the transverse roughmess 
loads for the smooth case 


load ratios 



CHAPTER 6 


EFFECTS OF CONSISTENCY VARIATION ON SQUEEZE Flli^ 

BETWEEN ROUGH SURFACES 

6.1 INTRODUCTION 

Sque©ze films are of practical interest in lubricaticai 
theory as they are representative of unsteady situations in 
hydrodynamic lubrication. The experiments conducted by 
Needs [ 1 ] and Fuks [ 2,3 ] to study the phencmenon of 
squeezing demonstrated a residual film after squeezing. 

Various explanations are offered for the creation of residual 
films. For example, Hayward and Isdale [ 4 ] suggested that 
apart from dirt and undissolved particles in the residual film, 
the surface asperities of the bounding surfaces could be offered 
as explanations- Derjaguin devised an experlmoital method to 
avoid any possible interference of dirt particles and observed 
an increase in viscosity when traces of polar or siirface 
active substances were present [ 5,6 ]. With the development 
of roughness models, squeeze films between rough surfaces 
have been studied [ 7,8 J. If certain rheological abnormalities 
Were obsered in thin films, a. natural question is that how 
these rheological anomalies affect the Rearing perfocmaace with 
surface asperities. Using micro polar fluid theory, Rrakash 
et.al [ 9 ] showed that the rheological anomalies like eihancem^t 

of viscosity in thin films ware rheological effects and could 
not be attributed to the presaice of asperities- In this Chapter, 
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we study the effects of surface roughness and rheological 
anomalies such as consistency variation across the film 
thickness by following a general model for consistency variation 
discussed in Chapter 5 together with Christenssi's theory 
of roughness in the case of power law lubricants. 


6.2 PARALLEL PLATES 

Consider the flow of a power law lubricant between 
two normally approaching infinite parallel plates with a 
normal velocity V (Fig. 6.l). The plates are of length 2d and 
are separated by a stochastic film thickness 2h.2h^ is the 
initial nominal film thickness of the systan just prior 
to squeezing . 

(i) ONE DIMENSIONAL LONGITUDINAL ROUGHNESS 


The stochastic form of Reynolds equation for this case 
of roughness can be obtained from eqn. (F>-13) and Fig. 6.1 
as 


d r _n 
dx^2n+l 




E(f ) (- - 
r m. 


■SE'P) i/m 

dx ' J ° ■ 


dh 


n 


dt 


( 6 . 1 ) 


where (f^) is given by 

(f c d-K-V") (h -a ) )/ ( 6 . 2 ) 

r n p 


where is apparent peripheral layer. 

To determine average pressure E(p) we use 


the following boundary 


conditions I 




\y 


Fig. 6*1 " Squeezing between two rough parallel ^plaits 
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dE(p) 

dx ^ ° at X = 0 
E(p) =0 at X = d 


(6*3) 


Integirating eqn. (6,1) and using conditions (6, 3), we obtain 
the average pressure E(p) as 


dh 


E(p) = m, (- (1 )q £.1 f - 

^ 1 n dt ' 'h/ ^E((f ))^ n+ 

1 r 


1+1 n+1 

-X 

n+1 


(6.4) 


The load W^(n,K,q,b) for the case of longitudinal roughness i 

I 

is given by I 


Wj^(n,K,q,b) = 2 / E(p) dx C6.5) 

o 

Using eqn. (6.4) , eqn. (6.5) becomes 

W^(n,K, q,b) = • 'A(_ -2n+l n)n ^^^q n+2 j 1 j 

L ' n+^ n dt h^ 


The squeeze time from the initial nominal film thickness 2hj^ to 
a subsequent nominal film thickness 2h2 is given by 

2m, 


. f „ , ^ / ^’^l )l/n f2n+l) (l_)q,/n in+2)/n ^1 J^^n 

t^Cn^K.q.b) .= (^+2 w:in,K7q,b)-^ ^ n ^ ^ ^ 


(6.7) 


The load components Wg (n,K^ q) and response time ■t^(n,K^q) 
for the smooth case are given by ( eqns, (3.16) and 

(3.17) ) 

/ V ^*^1 r 2n+l “^^n^n /l_)q n+2 _ . 1. „ 

Wg(n,K^q) = ^2 “n ’ 'dt '^i ^ (f^)^ 


( 6 . 8 ) 
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and 


tg(n,K,q) 


2m-[ 

= <TnV2rw,(„;K,iO->^^" 


dh 


/ — - ^ ^ 
h2 (f^)hi^2n+l-qJ7n 


(6.9) 


where (f^) is given by (eqn. (3_-ll)} 


-1/nx f , /, ^ (2n+l) /n 


(f ) = 1- {l~(l-a/h ) 
o n 


( 6 . 10 ) 


being the film thickness in the smooth case and a, peripheral layer. 


To study the effect of longitudinal roughness on load 
and response time we define the quantities W^g and t^g as 
follows I 


W. 


(n, K, q, b) -Wg (n, K, q) (F^)^ jj2n+l-q 


LB 


'LB 


Wg(n,K, q) 


{eCf^))}’^ 


-1 


tj^ (n,K, q, b) -tg (n,K, q) 


t_(n/K, q) 


1 dH 1 

= { / 

{eCf^)} 


dH 


" ^ 7 (F^)H^2n+l-qr/n 


1-1 


where 




(F ) = C d-K-V") (H -a^) + K-V" H<2"+lVnj ^q/n 

r n p 


(F 


) = 1_(1_K"^/’^) i l- (l-i/H^) 


(2n+l)/n 


/“} 


( 6 . 11 ) 


( 6 . 12 ) 


(6.13) 

(6.14) 


« = V«s' Hj = hj/hj^,ap = 


(6.15) 


= a/hj^, B = b/h^ 
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To study the effect of consistency variation on load and 
response time# we define the following quantities using 
eqn . (6 . 15 ) • 

W^(n,K,q,b) [e{ 

^LK = ^'(n, r;q;bT = [icF^'))]'' ' 

t, (n,K,q,b) 1 dH 1 dH 

"^LK ” ‘tj^(n/rrq/h^ ~ H e(P^) ' g(j^(2n+l-q)/nj 

(ii) ONE DIMENSIONAL TRANSVERSE ROUGHNESS 


{6.16) 

(6.17) 


In this case/ the Reynolds equation is given by 
(using eqn. (5.50) and Fig. 6.1)- 

f_[ 5,-P,, hf" fE(--" „ ^ j"""] = - 

dx^ 2n+l 1 /r: 


Following a procedure similar to that adopted in the case 
of longitudinal roughness/ we obtain the load and response 


(6,18) 


time as 


W^(n/ 


, 2n+l ;“'n)n (l_)q ^+2 j.f . 1. j (g.ig) 
,K,g,b) = ■„ at' 'h,' (f )“ 


ind ' / \ / 

2mj_ a/n,ln±lf(l ^n+2)/n 

t^(n/K,q/b) = (^2'"Wj:“n;K/q7by^ n H 


il^ 

^2 (fj’" 


( 6 . 20 ) 


urar ^ ^versB roughness and consistaacy 

To study the effects of transverse xuuy 

« -t-irnr- we define the following 
variation on loa,d and reppons®^ ' 


quantities 


= (E )nH2n+l-q- E( 1 -1 

TB o. n (F^)" ^ 




( 6 . 21 ) 


.]-l ( 6 .22) 



198 


E(1/(F^)’^) 

= / i E(1/(F^)’^ '^n/ ^ 

^2 ^2 
6.3 CIRCULAR PARALLEL PLATES 


(6.23) 

(6.24) 


The configuration of the. system is as shovni in 

Fig. 6.2. The plates are of radius R and separated by 

a film thickness 2h. 

(i) ONE DIMENSIONAL RADIAL ROUGHNESS 


The average Reynolds equation in this case can 


be obtained by using 


i — - ^ 

r dr'- ^n-tl 


h^’^ E((fp)r(- 


1 ^^^Pjl/nn 
^ ^ -1 


dh 

n 

dt 


/ 


I 


(6.25) 


The boundary conditions are given by 


0 at r = 0 (6.26) 

dr 

E(p) = 0 at r = R 

(a, 9 t;') nslra conditions given in eqn- (6.26) 
Integrating eqn* (6.25^ using 


we 


obtain the average pressure E(p) as 

E(p) = mp- 


(6.27) 


The Toaa W,(n,K,q,b) for the case of longitudinal t^igbness 
is given by 


w, (n,K,q,b) = s 2itr E(p)ar 

L O 


(6- 28) 




Fig. 6-2 Squeezing between two rough circular plates. 
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Using eqn- (6,27) in eqn. (6,28) we obtain 


n m. _ dh , , n 

W, (n,K,q,b) = (- 


The response time from an initial nominal film thiclcness 2h2^ to a 
subsequent nominal film thickness^ 2 h 2 is given by 


TC m,R' 


n+3 


tL(n,K,g,b) = (:;,A ..w..,S -E 


2n+l r , -)^’^-(l_)q/n.^l ( 6 - 30 ) 

W^(n,K,q,b) (n+3) ^1 ECf^)) 


The load Wg(n,K, q) and time of squeezing tg(n,K, q) for smooth 
case are given by 


( ^ 2n+l ^\>n /l_)qRn+3 1 

Wg(n,K,q) = ‘ ^'^3(- - 2n dt ^ \ 


and 


(6.31) 


^n+3 


ts(n,K„)= TrEi— 

S 2n w^(n,K,q)tY^+3) 1 ,^2 ^ 


(6.32) 


To study the effect of longitudinal roughness on load arid 
response time, we define the following quantities . | 


W 


LB 


= r(F )’^ H2n+l-q/ ( s((F ))’^)f ] -1 

Lon ^ 


1 dH, 


1 ]-l 

^LB = [ HTF-p / "" 


(6-33) 


• (6.34) 

\ 


(P ) and (F ) and non-dimensional schene 

where forms of function ( -^7 r 

are given by eqns. ( 6 . 13 ) - (6*i5) •• 

• oj. variation the following 

To study the effect of consistency 
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quantities are defined t 


«LK = ' 

( 2n+l— q) /n j 

/ E(P^)f 

(6.35) 

1 

dH 1 

dH 



eTf T / ^ 

r' H., 

n 

E(l^n+l-q)7m 

(6.36) 


(ii) ONE DIMENSIONAL CIRCUMPERENCIAL ROUGHNESS 

Stochastic form of Reynolds equation in this case can be 
obtained as 


1 d r n_ _ 

r drL2n+.l 


E(D- ,.)-Vn r { - r ^ " 

1 dr ^ -J dt 

^ (6.37) 


Integrating eqn. (6-37) using conditions (6.26) we can obtain 
E(p) . As in the case of longitudinal roughness, we can define 
similar quantities to study the effect of transverse roughness 
and consistency variation as follows I 

Wjg = (F^)" E £l/(F^)") -1 (6.38) 


'TB 


w, 


TK 


'TK 


[/ t E (1/(F^)")^-''" dH / )dHy ((P^) H^2n+l-g>/n)] _1 

"a (6.39) 

e!i/(f )" }/e i (6.40) 


/fE (1/(P )”))1/" dHy /fE (1/ dH„ 

2 ' Hj (6.41) 


6.4 ROLLER BEARING 

Consider squeezing between two rough identical rollers 
approaching each other with a normal velocity V. Ihe conf igur 
ation of the system is as shown in Fig. (6.3). 




Fig. 6*3 Squeezing between two rough 

rollers 
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(i) ONE DIMENSIONAL LONGITUDINAL ROUGHNESS 

The average Reynolds equation in this case is obtained 
fron eqn. (5.13) by putting U = 0 as 


dx"- 


n j^yn dE^)_)l/n-, _ (f: y,o) 

2n+l 1 r m, dx'^ -I 3t (6.42) 


where half of nomninal film thickness is given by 


^n = + X / (2R) 


(6.43) 

h^ being nominal inimmum j-ilin tnickness and_R the equivaJ-eht radius. 
Using the boundary conditions 


dE ( p) n 4- n 
”'dx ^ at X = 0 


E(p') 


(6.44) 

(6.45) 


0 at X = d 


where 2d is the length of the film strectch in the bearing 
clearance/ the average pressure E(p) can be determined as 


3 h . d n 

E(p) = -m, (- --£•)" (Ujl f *c 

1 n 3t h^ ^ {E((f^))} 

The load W, (n,K, q,b) is given by 


d 

W^(n,K,q,b) = 2 / E(p)dx 

“ o 

, Using eqn, (6.46), eqn. (6.47) beccmes 

dh *1 d ._n+l 

WL(n/K,q,b) = 


{ E ((fj)| 


dx 

n 


(6.46) 


(6.47) 


(6. 48) 
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The squeeze time tL(n,K,q,b) is given by 


i-J -r Of_// _ i 


’{ECf J)}"" 
r 


o 

(6.49) 


The load ratios W„ and t.„ are 


W. 


LB 


'j_^g oiiu. are given by 


[ ) ax/ 


O iE((F^))} 


n 


o (p )n g2n+l-q 

o n 


(6.50) 


where is given by eqn. (6.14), a = a and 

X = x/a , H^=hn/h^ = l+a^x2, H = h ,/h^,H^ ■ h^/h^. 


= a /(2Rh^), ip = Ep/hp, i = a/h^ 


(6.51) 


and 


'LB 


1 1 vri+ll 1 /v, i i V^^+i V 1 /r, 

I i r dn / f ( f ■— - ^ dH 

o . O' ^ jn H2n+l-q o 

2 O n 

(6.52) 


Hj ° { E((F •■))}'' 
^ r 


where 

X = x/b, H = h/h^ = H^+a^X^, 
a = u/b, , a = “ d / ( 2Rhj^ ) 

The quantities and tg|^,for the purpose of studying 

consistency variation of the lubricant are defined as 


W. 


LK 


1 ^+1 1 x^+1 

, -..p av ^j--(2n+l-qTA „n 


O i(E{F„)}' 


)r 


(6.53) 


(6-. 54) 


where the expressions (F^.), (-q^ non-dimensional scheme 

are given by eqns. (6.13 ), (6.14) and (6.5l) respectively 



205 


and 


r _ r s r •, 1/n ^ ^ 

^LK ” ^ -v.™. ... } . dH / / { / - -- , 


^2 ®{E((P^))}^ 


i /n. 

H 2 ' o { :< H ^ 2 n+l-q) /njjn 

(6.55) 


where non-dimensional schane is given by eqn. (6.53). 

(ii)ONE DIMENSIONAL TRANSVERSE ROUGHNESS 

The stochastic form of Reynolds equation for this 
type of roughness is given by 


a_r i h?/" { E(-i-. ))-l/n (. 1 a|(p),l/n., _ _ 

■I. Or, 4.1 ± (fp mdx j 9t 


(6.56} 


The load ratio and response time ratio t^^ are obtained 


TB 


as. 


w, 


TB 


= [ / E { '~^} dX/ / {(P } dX]-l 

O (F„) o " ° 

^ (6.57) 


where the expressions (P^), (P^) and relevent non-dimeisional 

quantities defined in them are given as in the case of 

and 


1 1 


^TB = ■»„/){/''' -X^n, n+l-q 


Hj.o 


<V> 


H O 
2 


(f^)“h; 
o n 


the expression (Pj_)r^E^) and non-dimensional scheme are 
given as in the case of tj^g* 


(6.58) 
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To study the effect of consistaicy variation on load and 
response time we use the following quantities 

B i- -^-^}dX/ / E{~ p- ^ • } dX 

o (F^)^. o „2n+l-q 


(6.59) 


^TK ~ i ^ ^ 

K_ r> /rn ’kB ^ O' 


H2 ° 


(F^P 


1 1 

f i f X^+1 E( Vn 


^2 ° 


j2n+l 


(6.60) 


where the function (f^), (f^) and non-dimensional scheme 

for and are given by those of and tj^ respectively. 


6.5 RESULTS AND DISCUSSION 

To study the effect of K on load ratio and in the 
case of parallel plates, ecps. (6.16) and (6.23) are evaluated 
for various values of n and depicted together in Fig. 6.4. it 
is observed that as in the case of squeezing for smooth bearing, 
an increase in K results in an increase in load capacity for all 
n for both the cases of longitudinal and transverse roughness ^ 
pseudoplastic lubricants are much more influenced by the consistency 
variation across the film thickness than the Newtonian and 
dilatant lubricants. The response time increases both for 
longitudinal and transverse roughness for increasing K as is 
evident from Pig, 6,5 (using eqns,6.55 and 6.60) . 

Fig. 6.6 gives the picture of response 
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time vs. roughness parameter B ieqcis. 6,5Z and 6.58). As 
B increases, the response time decreases for the longitudinal 
case and increase for transverse case. Further, the effect 
of roughness is to increas®the response time for increasing 
flow behaviour index. 

The effect of the roughness parameter on load ratio 
^LB ^TB presented in Table 1. 

Table 1 

Variation of load ratios and with roughness parameter B 



\b 


W 

TB 



B m=0 . 5 

n = 1.0 

n = 1.5 

n = 0.5 

n = 1.0 

n = 1. 

.1 -0.534 

-0.721 

-0.83 2 

0.005 

0.013 

0.024 

.2 -0.536 

-0.722 

-0.833 

0.024 

0,057 

0,107 

.3 -0.540 

-0.725 

-0.835 

0.057 

0.141 

0.274 

K=1.5, a 

= 0,6, q 

= 0.5 





As the roughness parameu^er B increases, the load 
increases in the transverse case and decreases in the 
longitudinal case for all n. The effect of surface roughness 
for fitted K on loa.d is more for dilatants compared to 
Newtonian and pseudoplastic lubricants. 

In the case of circular parallel plates, the load and 
res non se time ratios obtained for radial and circumferential 



208 


ona dimensional roughness are the same as those the 
longitudinal and transverse roughness for parallel plates. 
Hence, similar interpretations can be accorded to circular 
parallel plates. 

In the case of rollers, eqns. (6.54), (6.55), (6.59) 

and (6»60) are evaluated to determine the effect of K 
on load ratio and response time and depicted in Figs. 6.7 
and 6.8 respectively. Similar to squeezing for smooth 
case , the load and response time increase as K increases 
in the case of rough rollers- These bearing characteristics 
follow the same trend with regard to flow behaviour index 
as in the case of smooth rollers under squeezing. 


The important aspect studied in this Chapter is the 
effect of roughness in the lubrication of fluid with 
layers of different consistencies. The effect of roughness 
on load ratios and can be studied by evaluating 

eqns. (6.50) and (6.57). This is depicted in Fig. 6- 9. 

The affect of roughness is to decrease the load capacity 
for longitudinal case and increase it for transverse 
case for all n. As n increases, the load ratio increases 
in the transverse case, implying that the dilatants are 
capable of supporting more load roughness . The dilatant 
lubricants show an opposite trend with regard to load 
capacity for the longitudinal case- 
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Similar behaviour of the dilatants with regard to 
response time is observed in Pig. 6.10 using eqns- (6.62) 
and (6.58). 


Table 2 gives the picture of tine effect of thermal factor 
q on load ratios and and response time ratios tj^^ 

and t^j^. It is observed that for K > 1 as q increases the 
load ratio decreases for the longitudinal roughness whereas 
it increases for transverse roughness. However, the increase 
is very small. The response time ratio decreases as q increases 
for both the cases of roughness structures. 


Table 2 

Effect of thermal factor q on load ratios and 


and 

response 

time ratios t^^^^ 

and t^j^ for 

q 

.^^LK 


"‘^TK 

.1 

1.3736 ' 

1.3646 2.129 

2.113 

.5 

1.3721 

1.3651 2.114 

2.101 

.7 

1.3713 

1.3653 2.106 

2,014 

.9 

1.3705 

1.3656 2.09o 

2 .087 

K = 

1.5, :.n 

= 0.5, ip = 0.4, 

w 

II 

o 

• 

K) 


6.6 G0^7CLUSI0N 

In the presence of asperities, the effect of consistency 
variation across the film thickness for squeezing for parallel 
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infinite plates, circular plates and rollers is the same 
as that of the smooth case for the respective configurations. 
The effect of roughness with specified consistency .variation 
parameters is more felt in the case of dilatants compared to 
Newtonian nnd pseudoplastic fluids with regard to load 
capacity and response time- The effect of the thermal 
factor parameter is to decrease the load ratio for longitudinal 
roughness and to Increase the ratio for transverse roughness 
when the peripheral lay-er has greater consistency than that 
of the central layer. However, the increase due to the effect 
of q is very small. 




> K 


The effect of consistency ratio K on load 
ratios and Wjj^for different values o1 
flow behaviour index n for parallel plate 



'LK/TK 



0*1 i-. X _. . ■ ,w l W _v- . K r :., I- . 

0*4 0-7 10 1-3 1-6 

K 


Frg. 6*5 Effect of consistency ratio K on respone 
t ime ratios t^j^ tfK for various values 
^ index n for pdraUei plates 





Fig.S-6 Effect of roughness B on the response 

time ratio t^Q and tjB for various values 
of flow behaviour index n for paraUelplates 






Fig. 6*7 Effect of consistency ratio K on load 
ratios and for various values 
of flow behaviour index n for rollers 





Fig. 6*8 Effect of consistency ratio K on 
response time ratios tj_|^ and 
for various values of flow behaviour 
index n for rollers 



Fig. 6-9 Effect of roughness parameter Bon 

load ratios W^^gand W|g for various values 
of flow behaviour index for rollers 
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a 


P 


b 


e( ) 
h 


^ 1 ' ^2 


n 


Hi,H2 




d 

d 




n 

Vr P]_/ P2 

q 

r 

R 


t 


apparent peripheral layer film thickness 
measured from average film thickness 
non-dimensional quantity corresponding to 

^P 

roughness amplitude measured frc«n average film 
thickness 

expectation of ( ) 

stochastic film thickness 

initial and subsequent film thickness 

nominal film thickness 

roughness part of film thickness 

non-dimensional quantities corresponding 

to h^ and h2 respectively 

non-dimensional quantities corresponding to 

h^/hg , ^ respectively 

half the length of rectangular plate or 

half of film stretch in roller bearing 

consistency indices 

flow behaviour index 

hydrodynamic pressure 

thermal factor 

radial distance 

radius of the circular plate or equivalent 

radius of rollers 

time 
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t^p^response time ratios defined for longitudinal 
and transverse cases of roughness and for 
consistency variation across the film 
thickness for these cases. 

V normal velocity 

Wg load capacities for longitudinal / transverse 

roughnesses and smooth case 
^TK capacity ratio defined for longitudinal 

and transverse roughnesses and for consistency 
variation for these cases respectively 



CHAPTER-7 


EFFECTS OF CONSISTENCY VARIATION IN ELASTOHYmODYNAMIC 
LUBRICATION OF ROUGH ROLLING SURFACl^ 

7.1 introduction 

Heavily loaded machine ccmponents such as gears, 
cams, rollers etc. operate essentially in elastohydrodynamic 
regime. The assumption that bearing surfaces are aiKJoth 
is no longer valid whoi the bearing operating ccxiditions 
involve small film thicknesses as in. the case of eM 
lubrication. Under this cond-xt-ton- the surface asperities 
may interfere and be comparable with the film thickness 
of the operating lubricant and significantly affect the 
bearing performance. The topographical structure, 
deformation mechanism and participation trend Of these 
asperities in lubrication, particularly in the ehd regime 
constitute one of the major studies in Tribology- 
Fowles [ 1 ] studied the ehd lubrication betweei identical 
sliding asperities. Lee and Cheng [ 2 ] studied 
the effect of single asperity on the film and pressure 
distribution during its entrance into an ehd contact. 

Tallian [ 3 ] studied pressure and traction rippling in 
ehd contact of rough surfaces. 

Att 0 n:ipts were made to include the rheological properties 
of the lubricant to determine the characteristics of rough 
bearings. Furey [ 4 ] pointed out that lubricant viscosity 
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and roughness effects had a significant interaction in thin 
film lubrication. Using micropolar fluid theory, Prakash 
et.al [ 5 j concluded that the rheological anomalies cb>served 
in thin films could not be attributed to the presence of 
asperities of the bounding surface, Christensen [ 6 ] 
drew attention to the increase in temperature in fluid region 
possibly by the interaction of asperities in the inlet zone 
of ehd regime. But, Cheng and Dyson [ 7 ] pointed out that 
in the inlet zone of ehd regime ,the isothermal theory could 
be sufficient for the rollers operating with realistic values 
of coefficient of friction and upheld the validity of the assvjm- 
ption of equal tenperature for bearing surfaces and fluid 
region. These studieS/ in geieral, do not take into consideration 
the rheological anomalies such as enhanced viscosity in 
the convicinity of the bearing surface- 

In this Chapter, we analyse the el as to hydrodynamic 
lubrication of rollers by considering the surface roughness 
and consistency variation of the operating power law lubricant- 
Consistency variation is accounted through the model suggested, 
in Chapter 5 - The ehd analysis is limited to the inlet zone j 

i 

as the nominal film thickness, an aspect of primary interest i 
to the design engineer under, heavily loaded conditioi is j 
completely determined by performing such an aii.^ysis [ 8 ]. 
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7-8, ehd lubrication of rough rollers 

The problan considered is that of the ehd lubricaticai 

of rough rollers with a power law lubricant. The ccaisistency 

of the operating lubricant is considered to vary along 

as well as across the film thickness. The roughness effects 

are studied on the pressure distribution in the inlet zone 

following Christensen's approach [ 9 ]. To perform inlet 

analysis we resort to Grubin's theory [ 10 ]. Thus, the 

shape of the deformed surface at the contact is Hertzian, 

and in the main load carrying area the surfaces are separated 

by a parallel fluid film of mean nominal film thickness h . 

m 

It has been pointed out by Jackson and Caaeron [ 11 } and 
Kaneta and Cameron [ 12 ] that under pure rolling cCTiditions 
the surface asperities do not deform in the elastohydrodynamic 
lubricated (ehl) contact and Hertzian zone is nearly flat. 
Keeping this in view, in this Chapter/ the elastic deformation 
of the asperities just prior to Hertzian zone is assumed 
to be negligible [ 13 ] . Under typical conditions of ehd 
lubrication, contacting bodies are much more flexible than 
either the oil film or surface asperities and deformation 
of bearing surfaces can be calculated by the Hertzian 
theory [ 14 ]. Recently, Prakash and Czichos | 15 ] pointed 
out that in the absence of asperity interaction, it is 
sufficient to analysis the inlet region in the ehl contact- 
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The geometiry of the iubjclcated rollers deforming 
under the action of heavy load is depicted in Fig.7.1. 

In the ehd contact, the reduced pressure g* and the 
exponential variation of consistency with pressure are 
defined as £ 16 ] 

-K- ““Pp 

q = (l-e ®)/a (7.1) 

“Pe 

mi = ra^e (7.2) 

where a is the pressure-consistency coefficient and m 

o 

is the ambient value of the consistency mi, appearing in 
eqn. (2.14). mi can be treated to vary with pressure 
isothermally as the tanperature effects can be incorporated 
in the parameter q in the definition of m in eqn. ( 2.14). 

In fact, the general expression for consistency m can be 
given as [ 17 ] 

m = f(p)(h/hi)'^ (7>-3) 

(i) ONE DIMENSIONAL LONGITUDINAL ROUGHNESS 

Consider the deformation of rollers under heavily 
loaded condition when the peripheral layer thickness 
a is less than half of nominal film thickness h^^. The configurat- 
ion of deformation of rollas'ls depicted in Fig. 7.1, where 
2d is the width of Hertzian region. Pressure is maximum 
just outside the parallel width of the lubricant film, i.e. \ 
at X = -d. \ 




Fig.T'l Elastohydrodynamic lubrication between 
two rough rollers considering consistency 
variation hf^^^Zap 
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lAJith the help of eqns. (7.1). t7.2) and ('5.1'^), the modified 
Reynolds eguation for the case of longitudinal roughness can be 
written as 


d r 

dx ^ ^n+r 


ih hf E (f ) (i» ^ 


-oo < X < — d 


(7.4) 


where h^^ denotes the nominal inlet film thickness which is 
sufficiently away fran the contact zone where q is measured * d 
is given by eqn. (2.50) ^ 

^ = 2[Wg Rg^/E']. (7.5) j 

where ^ /2)/E' = il-v^}/E^, R^ = r^/2‘ (7.6) 


Wg /Rg , E^ and v being load in the ehd lubrication, radius 
of the equivalent roller. E^ Young's modulus and v Poisson's ratio 
of the rollers. 

The expression (f^) is given by (eqn:. 5.10) 

(f^) ={ (1-K~^/’^) (h^-2ap) + K“^/n j^(2n+l)/nj /^q/n 

The smooth part of the film thiclcness h outside the band of 

n 

contact width is given by (using eqn 42 .S 4 ) 


h^ = h^ + d^/(2R)[ Ix/dl (x^ ) ~l-ln((| xl d ) -s- ^ (x^/d^-D] (7.8) 

where h - h^^ is the minimum nominal film thickness 
The boundary conditions to determine pressure E(q'^^) are the 
following 1 


E (q '') = 0 at x = - 00 (7.9) 

Mq::l= 0 at X = -d (7.10) 

dx 


Integrating eqn. (7.4) with conditions (7,9) and (7.10) . we ' 
obtain average pressure E(q*) for the inlet zone. Denoting E(q'’'‘) 
at X = -d for the longitudinal roughness case as q* we have 
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q» = [E(q*)] 


(7,11) 


x=~d 

= (2^:Ll)n m 

2n ■ ^ V ^ ... dx 

^ iEKf^))}^ 

The ireduced pressure q obtained for the smooth case in Chapter 2 

(eqn. 2,59) at x = -d, denoted by is the giv® by (eqn. (2,59)} 

<Is = [<3 ]x=-d 

„ 1 „ -<2 <h_-h_)" 

dx (7.12) 


= (?ntl)n jSn+l m u" ( 1)<3 "/ 

1-00 .n 2n+l-q 

o n 


where 


(f„) = l-Cl-K-V") {l-(l_ 2 a/h„)' 2 n+i>/n , 


-1/n, 

f 1 ( 1 O /"U ^ X / /n ,1 ^ 

To study the effect of longitudinal roughness on the pressure 
distribution in the inlet region, we define - a quantity given by 


<• 


-{. ^^"^5 “1 -H )’^ 

o = _ I f ^-.J^ / r n m jv , T 

L -CO '■ (f )’^H2^+l-g 


{ E((F^)) ]' 


n 


(7.14) 


where 


H = (hE')/(2Wg) = H^+H^ , = (h^E^ )/(2Wg),H^ = (hgE')/(2Wg) 


(7.15) 


= (h^EO /{2VJg), i = (aE') /(2Wg), = (a^^'’ ) / (2Wg) , B =(bEM/(2Vi 




^'2 


= H^, +i X\ X^-l-ln(IXl ’ X"^-1) (7. IS) 

(F ) = 1-(1-K“^^’^) { (l-2i/H ;|^2n-!-l) /n (7.17) 

(Fj.) = { (l-K-^'")(H^-2i )^2n+l)/n^^-l/n jj(2n+l)/n, ^j^/n 

The result obtained in eqn. (7.13) for q£ is valid for h^^ 
To obtain the factor for 'h^2ai^ , we divide the inlet 


zone' into two regions -^ < x < -dj^ where h^ ^ 2a, 


o 


and 
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-di < x<^d, where \<2a.^ and = 2ap= hg, at x = -d^ 

(Fig, 7,2). 

The governing differential ec^i. for reduced pressure = q* 
in the region -00 < jc < -d^ can be written from eqn. (7,4) 
as 


^ r_2n /I) (2n-i-l)/n ,g/n 

dx L 2n+l4^ dx ^ J = 


dh 


u 


n 

dx 

(7.19) 


where represents the reduced pressure q for h > 2a 
In the region -d^^ i ^ i the consistency of the film 
is Km (h/hj^ ) '^ over the entire region and we can write 
the differential eqn. for pressure in this region as 


fe-in+j‘l> g(j^2n+l-g)/nj 




dE(q*) 

... 2_3l/n 

Km dx ‘ 
n 


= U dx 


(7 ,20-) 

where q 2 denotes the reduced pressure q"^ for hn< ^a^* 

To determine reduced pressure from eqns. (7.18) and (7,19) , 
we use the following conditions 


E(q^) =0 at X 

dE(q*) _dE(q|) 

~dx 'djT 


-CO 


at X = -d^ 


dE (<52^ 


dx 


= 0 


at x = -d . ^ 


;7 .21 ) 
K 7.221 

(7.23) 


Integrating eqn. (7 .19) we get 




(hn+h^) 
i E((f jJ)} 


n 


n 


(7.24) 




> X 


Fig. 7*2 F lastohydrodynmic lubrication between two 
rough rollers considering consistency 
vGtiation hm^20p 



231 


where is the constant of integration to be determined 

subject to the condition h. = 2a at x - -d. 

n 

Integrating eqn. (7.19) subject to condition (7,23) we 
have 


( ''i XX 

(-2p+l)n 22n+l ^ j 

2n o h^ |E(j^(2n+l-q)/n)p 


where hn = hi at x = -d. 
in 


Using condition (7.21) in eqns. (7. la) and (7.19), h. can 


be determined as 




E(lE )) 
r 


^ntl-q) /n^ h^=2; 


(7. 26) 


n p 


As in the case of 1r^> 2a , we obtain the reduced pressure 
^ ^ ) at X = -d as 


4 = 

-d 

(2n+l^n ^2n+l ,,n,l vq. . ^ ^ 

= (“-or") 2 ru u (^) ^[ I _ - - dxt 


° ^1 -CO iE((f^))}'^ 


d (h -h ) 

K / - . ; — dx 


(7.27) ^ 

i 

Using the dimensionless quantities given in eqn. 7.15 I 

^ _ • 

we define the pressure ratio factor Qt for Hn < 2a as 

b — p 




1 (H^+Kr (H -H D -H ) 

— - — dX +K •!1 '-r-n .. Tn _ dx} / I ^ -'6 dx 

D Mf^))}^ “°1 £ E((P^))i^ 


T (K -H 

K ^ ^ . dx }] 

/ j^2n+l-q -* 

_D^ “n 


(7.28) 



where and 


H. = 


= = KVn(25 H ) [- 


e((f_)) 


J-T AAt/AX 

L 1/ m P ^ 


E(H< 2 n+l-q)/n) H^=25p 


.2i 


(ii) ONE DIMENSIONAL TRANSVERSE ROUGHNESS 


(7.29) 


^or this type of roughness, the differential equation 
to determine reduced pressure for the case hj^^^ 2a can be 
rewritten from eqn. (5.49) as 

^r.J^„»(i.)^2n+l)/n-,q/nj , 1 ) . -1/n A dE(g*) 1/n n 
dxL 2n+l4^ ^1 jn/^ 'dx J 

dhj^ 

^ ^ "dx ^ ™ 

Using the boundary conditions (7. 21)- (7.23 ) and following 
a similar procedure adopted in the case of longitudinal 
roughness, we obtain the reduced pressure q* at x = -d 
as 


q 


* 


[E(q*)] 


x=-d 


^(2n+l)n 22n+l^ ^n^l (h^-h^)’" E(--1- - .dx (7.31 ) 

n o Jri-i A 

1 -oo r 


and pressure ratio factor as 

-1 . ^ -1 (H -H )"■ 

_ r f ^n.,./ 1 ^ n m 


q: = [ 


;) dX/ 


^ ■- -co-n-ta- J^\^Jn^2n+l-q 


■■■■ dX]-l 


(H: -H ) e( 

-oo - ^^ 3 ^) 

In the case of h^'lSa, as in the longitudinal case, the inlet 
region is divided into two region where h > 2a and '!^i£2a and 


n- ^ n 

the set of equations to determine reduced pressure in these regions 
can be obtained as 
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= U 


.‘^'v 

dx 


(fp' 


m dx 
o 


^ yJ. /11-j 


-OO < X < -d. 


r_2n /l^ (2n+l)/n , q/n , 1 _s,-l/n , 1 

&'’2n+l^r 


1 

dh^ 
= U-d-x 


(7.33’) 




dx 


(7.34.) 


where and q* represent reduced pressure q"^ for the 
regions defined in eqns. (7,3 3) and (7.34) respectively. 
Using the boundary conditions (7. 21 ) - (7.23 ) we can obtain 
the reduced pressure at x = ~d as 

^(2p+ljn22n+l^y(l / (h^+hj,)"E( U dx 

1 -00 vf ) 




(7.3 5) 


where 


eC 


) 


hj = Kl/"(2ap-ht„)[ 


jj2n+l-q 

E (-• : 

(f‘3,)’^ 


■t/'rt 

~] -2 a 

h. =2a„ ) 


n . r- 

and pressure ratio factor Q as 


C7 •3 6) 
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-D, 




-oo 


n ~ -D. (H -H 

1 „ -1 1 n m 

n j^2n+l-q J- )n^2n+l-q 


1 Xh -a 

K / --dX } -1 

■’ 2n+l-q 

“^1 


(7.37) 


H is given by 
T 


g(l/H2n+l-q) _ 


H = -H )[ - — n”] ' 

T Pm'- e(1/((F H = 25 P 

^ r n p ■ 


il/n 


(7.33) 


7.3 RESULTS AND DISCUSSION 

The. pressycf® ratio factors Qj^ and can be evaluated 
from eqns. (7.14) and (7,32) respectively in the case of 
^ 2ap for various values of roughness parameters,^ the 

consistency ratio K and thermal factor q. Fig. 7*3 represents 

* * 

the variation of and xvith B for various values of n. 

It is observed that as B increases, increases ^ the increase 
is enhanced as flow behaviour index increases. With increasing 
value of B^ ■ decreases and the decrease is ehanced with 
increasing n. However, the decrease is not significant. 

The results are in good qualitative agreement with those given 
in the reference [ 30]. 

« * 

Fig. 7,4 gives the picture of the variation of and 
with K in the presence of asperities. For very high and law 
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values of K, the dilatants are much more influenced than 
Newtonian and pseudoplastic lubricants. The influence results 
in increasing inlet pressure for K > 1 and this pressure 
decreases for K < 1 for dilatants. However, the variation 
is small, indicating thereby the rheological effects cannot 
be attributed to the presence of surface asperities. This 
is iTii- conformity with the results of Prakash et. al [ 5 ] 

who studied slider bearing using a micropolar fluid. 

* * 

Vilhen < 2ap , the pressure ratio factors and are 
evaluated by using eqns . (7.28) and (7.37) and depicted in 

Figs. 7.5 and 7.6. In the pr essence of surface asperities, 
the effect of K is to increase the pressure ratio factor for 
both longitudinal and transverse cases for all K (Fig. 7*5). For 
K > 1, the increase of this factor due to the effect of K 
is more for dilatants for both the roughness patterns. For 
K< 1, the trtsnd is reversed. The effect of the roughness 
paraincter B is to increase the pressure ratio factor for 
transverse roughness and decrease it for longitudinal roughness 
(Fig. 7. 6) . 

These results may be interpreted in terms of nominal film 
thickness at the inlet. An increase or decrease of pressure 
ratio factor will result in an increase or decrease in the 
nominal inlet film thickness. Thus, it can be concluded that 
the transverse roughness increases the nominal inlet film 
thiclaai;;ss anri the longitudinal roughness decreases it. Larger 
the value of flow behaviour index n, greater the rate of incroas, 
or dccronSG of nominal inlet film thickness. 
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7.4 CONCLUSION 

In the ahd lubrication of rough rollers, an increase in 
the consistency in the peripheral layer results in an increase 
in the pressure ratio factor or inlet film thickness. In case 
of transverse roughness, the pressure ratio factor increases 
as the value of the roughness parameter increases, for both 
the cases when the total peripheral layer thickness on the 
two rollers is greater or less than the minimum nominal 
fili.i thickness i In the presence of consistency variation 
across the film thickness it can be concluded that the dilatants 
are more effective than pseudoplastics in the case of transverse 
roughness as compared to the case of longitudinal roughness. 

The effect of the thermal factor is to the decrease the pressure 
ratio fvactor for transverse case and to increase it for 
longitudinal case. It is concluded that the rheological 
anomalies in the lubrication process cannot be attributed to 
the presence of surface asperities for non-Newtonian lubricants 
characterized by the pov;er law model. 
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NOMENCLATURE 


a 

a 

d 

E ( ) 
h 
h 


^1 


peripheral layer thickness 

non-dimensional quantity corresponding to a 
roughness amplitude measured from the average 
film thickness 
expectation of ( ) 

film thickness 
inlet film thickness 
minimum film thickness 
nominal film thickness 

roughness part of stochastic film thickness 
H, Hn, Hs non-dimensional quantities corresponding to 

h/hj^,h^,h^ and h^ resply 
consistency ratio 
half of hertizian width 


h 


h 


n 


h 


K 

d 


% 

D 


m/m^ym^ 

n 

Pe 

-){■ -ss- 

q /qi^q2 

O''^ 0*'^ ' 

^ *"r 


the point at which h^ =2a 
ratio = d/d^ 
consistency indices 
flow behaviour index 
elastohydrodynamic pressure 
reduced pressures 
pressure ratio factors 



